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Abstract 



' We formulate the basic postulate of pre-big bang cosmology as one of "asymp- 

o 

00 

' perturbative solution of the tree-level low-energy effective action. Such a 



totic past triviality", by which we mean that the initial state is a generic 



past-trivial "string vacuum" is made of an arbitrary ensemble of incoming 
gravitational and dilatonic waves, and is generically prone to gravitational 



instability, leading to the possible formation of many black holes hiding sin- 
gular space-like hypersurfaces. Each such singular space-like hypersurface of 



^ . gravitational collapse becomes, in the string-frame metric, the usual big-bang 

t = hypersurface, i.e. the place of birth of a baby Friedmann universe after 
a period of dilaton-driven inflation. Specializing to the spherically-symmetric 
case, we review and reinterpret previous work on the subject, and propose a 
simple, scale- invariant criterion for collapse/inflation in terms of asymptotic 
data at past null infinity. Those data should determine whether, when, and 
where collapse/inflation occurs, and, when it does, fix its characteristics, in- 
cluding anisotropics on the big bang hypersurface whose imprint could have 
survived till now. Using Bayesian probability concepts, we finally attempt 
to answer some fine-tuning objections recently moved to the pre-big bang 
scenario. 
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I. INTRODUCTION AND GENERAL OVERVIEW 



Superstring theory (see []T| for a review) is the only presently known framework in which 
gravity can be consistently quantized, at least perturbatively. The well-known difficulties 
met in trying to quantize General Relativity (GR) -or its supersymmetric extensions- are 
avoided, in string theory, by the presence of a fundamental quantum of length ~ (a')^/^. 
Thus, at distances shorter than Is-, string gravity is expected to be drastically different from 
-and in particular to be much "softer" than- General Relativity. 

However, as was noticed since the early days of string theory 0, a conspicuous difference 
between string and Einstein gravity persists even at low energies (large distances). Indeed, 
a striking prediction of string theory is that its "gravitational sector" is richer than that of 
GR: in particular, all versions of string theory predict the existence of a scalar partner of the 
spin-two graviton, i.e. of the metric tensor g^^, the so-called dilaton, 0. This field plays a 
central role within string theory |l|] since its present vacuum expectation value (VEV) 0now 
fixes the string coupling constant g"^ = e"^""", and, through it, the present values of gauge and 
gravitational couplings. In particular, it fixes the ratio of ig to the Planck length £p ~ G^/^ 
by ip ~ gig- The relation is such that gauge and gravitational couplings unify at the string 
energy scale, i.e. at E ~ Mg ~ gMp ~ (a')"^/^ ~ 3 x 10^'' GeV. It thus seems that the 
string way to quantizing gravity forces this new particle/field upon us. 

We believe that the dilaton represents an interesting prediction (an opportunity rather 
than a nuisance) whose possible existence should be taken seriously, and whose observable 
consequences should be carefully studied. Of course, tests of GR 0] put severe constraints on 
what the dilaton can do today. The simplest way to recover GR at late times is to assume 
^ that gets a mass from supersymmetry-breaking non-perturbative effects. Another 
possibility might be to use the string-loop modifications of the dilaton couplings for driving 
(j) toward a special value where it decouples from matter 0. These alternatives do not rule 
out the possibility that the dilaton may have had an important role in the previous history 
of the universe. Early cosmology stands out as a particularly interesting arena where to 
study both the dynamical effects of the dilaton and those associated with the existence of a 
fundamental length in string theory. 

In a series of previous papers 0, a model of early string cosmology, in which the 
dilaton plays a key dynamical role, was introduced and developed: the so-called pre-big 
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bang (PBB) scenario. One of the key ideas of this scenario is to use the kinetic energy of the 
dilaton to drive a period of inflation of the universe. The motivation is that the presence of 
a (tree-level coupled) dilaton essentially destroys the usual inflationary mechanism [|T^: 
instead of driving an exponential inflationary expansion, a (nearly) constant vacuum energy 
drives the string coupling g = e*^/^ towards small values, thereby causing the universe to 
expand only as a small power of time. If one takes seriously the existence of the dilaton, 
the PBB idea of a dilaton-driven inflation offers itself as one of the very few natural ways 
of reconciling string theory and inflation. Actually, the existence of inflationary solution in 
string cosmology is a consequence of its (T-) duality symmetries . 



This paper develops further the PBB scenario by presenting a very general class of pos- 
sible initial states for string cosmology, and by describing their subsequent evolution, via 
gravitational instability, into a multi-universe comprising (hopefully) sub-universes looking 



like ours. This picture generalizes, and makes more concrete, recent work [0, about 
inhomogeneous versions of pre-Big Bang cosmology. 

Let us first recall that an inflation driven by the kinetic energy of forces both the 
coupling and the curvature to grow during inflation 0. This implies that the initial state 
must be very perturbative in two respects: i) it must have very small initial curvatures 
(derivatives) in string units and, ii) it must exhibit a tiny initial coupling g^ = e*^*/^ <^ 1. 
As the string coupling g"^ measures the strength of quantum corrections (i.e. plays the role 
of h in dividing the lagrangian: C = g~^ C), quantum (string-loop) corrections are initially 
negligible. Because of i), a' corrections can also be neglected. In conclusion, dilaton-driven 
inflation must start from a regime in which the tree-level low-energy approximation to string 
theory is extremely accurate, something we may call an asymptotically trivial state. 

In the present paper, we consider a very general class of such "past-trivial" states. Actu- 
ally, perturbative string theory is well-defined only when one considers such classical states 
as background, or "vacuum", configurations. For the sake of simplicity the set of string 
vacua that we consider are already compactified to four dimensions and are truncated to 
the gravi-dilaton sector (antisymmetric tensor and moduli being set to zero). Within these 
limitations, the set of all perturbative string vacua coincides with the generic solutions of 



the tree-level low-energy effective action ||n 



Ss = — fd^xVG e-'^lRiG) + Gf"'d.(l) d^<l)] , (1.1) 
a' J 
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where we have denoted by G^^^ the string-frame (a-model) metric. The generic solution is 
parametrized by 6 functions of three variables. These functions can be thought of classically 
as describing the two helicity— 2 modes of gravitational waves, plus the helicity— mode of 
dilatonic waves. [Each mode being described by two real functions corresponding, e.g., to 
the Cauchy data (0, 0) at some initial time.] The same counting of the degrees of freedom 
in selecting string vacua can be obtained by considering all the marginal operators (i.e. 
all conformal-invariance-preserving continuous deformations) of the conformal field theory 
defining the quantized string in trivial space-time. We therefore envisage, as initial state, the 
most general past-trivial classical solution of (|1.1|) , i.e. an arbitrary ensemble of incoming 
gravitational and dilatonic waves. 

Our aim will be to show how such a stochastic bath of classical incoming waves (devoid 
of any ordinary matter) can evolve into our rich, complex, expanding universe. The basic 
mechanism we consider for turning such a trivial, inhomogeneous and anisotropic, initial state 
into a Friedmann-like cosmological universe is gravitational instability (and quantum particle 
production as far as heating up the universe is concerned f^). We find that, when the initial 
waves satisfy a certain (dimensionless) strength criterion, they collapse (when viewed in the 
Einstein conformal frame) under their own weight. When viewed in the (physically most 
appropriate) string conformal frame, each gravitational collapse leads to the local birth of a 
baby inflationary universe blistering off the initial vacuum. Thanks to the peculiar properties 
of dilaton-driven inflation (i.e. the peculiar properties of collapse with the equation of state 
p = e), each baby universe is found to contain a large homogeneous patch of expanding 
space which might constitute the beginning of a local Big Bang. We then expect each of 
these ballooning patches of space to evolve into a quasi-closed Friedmann universe Q. This 
picture is sketched in Fig. H and Fig. ^ In order to study in detail this scenario, we focus, 
in this paper, on the technically simplest case containing nontrivial incoming waves able 
to exhibit gravitational instability: spherically symmetric dilatonic waves. This simple toy 
model seems to contain many of the key physical features we wish to study. 

Before entering the technicalities of this model, let us clarify a few general methodolog- 



^Our picture of baby universes created by gravitational collapse is reminiscent of earlier proposals 
|15|, M], 1 17], but differs from them by our crucial use of string-theory motivated ideas. 



TIME 

FIG. 1. Symbolic sketch of the birth of many pre-big bang bubbles from the gravitational insta- 
bility of a generic string vacuum made of a stochastic bath of classical incoming gravitational and 
dilatonic waves. Each local Einstein-frame collapse of sufficiently strong waves forms a cosmologi- 
cal-like space-like singularity hidden behind a black hole. The parts of those classical singularities 
where the string coupling grows inflate, when viewed in the string frame, and generate ballooning 
patches of space (here schematized as the stretching of one spatial dimension ) which are expected 
to evolve into many separate quasi-closed Friedmann hot universes. 

ical issues. One of the goals of theoretical cosmology is to "explain" the surprisingly rich 
and special structure of our universe. However, the concept of "explanation" is necessarily 
dependent on one's prejudices and taste. We wish only to show here how, modulo an "exit" 
assumption [^, one can "explain" the appearance of a hot, expanding homogeneous universe 
starting from a generic classical inhomogeneous vacuum of string theory. We do not claim 
that this scenario solves all the cosmological "problems" that one may wish to address (e.g., 
we leave untouched here the monopole and gravitini problems). We do not try either to 
"explain" the appearance of our universe as a quantum fluctuation out of "nothing". We 
content ourselves by assuming, as is standard in perturbative string theory, the existence of 
a classical vacuum and showing how gravitational instability can then generate some inter- 
esting qualitatively new structures akin to those of our universe. In particular, we find it 
appealing to "understand" the striking existence of a preferred rest frame in our universe, 
starting from a stochastic bath of waves propagating with the velocity of light and thereby 
exhibiting no clear split between space and time. We shall also address in this paper the 



question of the naturalness of our scenario. Recent papers |T^, [|T9l have insisted on the 
presence of two large numbers among the parameters defining the PBB scenario. Our an- 
swer to this issue (which is, anayway, an issue of taste and not a scientifically well posed 
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FIG. 2. A representation of pre-big bang bubbles similar to that of Fig. but in 2+1 di- 
mensions. The different horizontal planes represent different instants in the evolution from the 
asymptotic trivial past to the Friedmann phase. Two inflationary bubbles characterized by two 
different initial horizon sizes (both large in string units) are shown to lead to Universes of very 
different homogeneity scale at the time at which the Hubble radius reaches string-scale values 
(Is = 0(10-32 cm)). 
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problem) is twofold: on the one hand, we point out that the two large numbers in question 
are classically undefined and therefore irrelevant when discussing the naturalness of a classi- 
cal vacuum state; on the other hand, we show that the selection effects associated to asking 
any such "fine-tuning" question can render natural the presence of these very large numbers. 

II. ASYMPTOTIC NULL DATA FOR PAST-TRIVIAL STRING VACUA 

Motivated by the pre-Big Bang idea of an initial weak-coupling, low-curvature state we 
consider the tree-level (order g'"^) effective action for the gravitational sector of critical su- 
perstring theory, taken at lowest order in a', Eq. ( |1 . 1| ) . We set to zero the antisymmetric 
tensor field B^^, and work directly in 4 dimensions (i.e., we assume that the gravitational 
moduli describing the 6 "internal dimensions" are frozen). Though the physical interpre- 
tation of our work is best made in terms of the original string (or cr-model) metric G^jy 
appearing in Eq. ( |1 . 1| ) , it will be technically convenient to work with the conformally related 
Einstein metric 

g^^ = e-('^-<^-w) Q^^ ^ IQ^Q ^ ^/ g<^now _ (^2.1) 

In the following, we will set G = 1. In terms of the Einstein metric Qfj^iy, the low-energy 
tree-level string effective action (|1.1|) reads 



(2.2) 



The corresponding classical field equations are 



R^u = ld^(l)d,<P, (2.3) 



□ = V^V^0 = O, (2.4) 



with Eq. (2^) actually following from Eq. ( p.3|) thanks to the Bianchi identity. As explained 



in the Introduction, we consider a generic solution of these classical field equations admit- 
ting an asymptotically trivial past. Such an asymptotic "incoming" classical state should 



^We will use the signature (— , +, +, +) and the conventions R^^p^ = dpT^^ — ■ ■ ■ , R^iu = R'\pi,- 



7 



allow a description in terms of a superposition of ingoing gravitational and dilatonic waves. 
The work of Bondi, Sachs, Penrose and many others in classical gravitation theory indi- 
cates that this incoming state should exhibit a regular past null infinity X~ , and should be 
parametrizable by some asymptotic null data (i.e. conformally renormalized data on X~). In 
plain terms, this means the existence of suitable asymptotically Minkowskian ^ coordinates 
(x^) = {t, X, y, z) (which can then be used in the standard way to define polar coordinates 
r,9,Lp and the advanced time v = t + r = t + (x^ + + z'^Y^'^) such that the following 
expansions hold when r — > oo, at fixed v, 6 and y?: 

^(,^)=^„ + Z(!^ + „fiV (2.5) 



= + (2.6) 



The null wave data on 1~ are: the asymptotic dilatonic wave form f{v,6,(p), and the two 
polarization components f+{v,6,ip), fx{v,6,(p), of the asymptotic gravitational wave form 
ffii^iv, 0, ^p) 0. [Introducing a local orthonormal frame e(i), e(2) on the sphere at infinity one 
usually defines /+ = |(/(i)(i) — /(2){2)), fx = /(i){2)-] The other o(l/r) pieces in the metric 
are gauge dependent, except for the "Bondi mass aspect" moo{v,9,ip) (defined below in a 
simple case) whose advanced-time derivative is related to the (direction-dependent) incoming 
energy fluxes: 

|- m^{v, e,if) = \ {d, ff + \ {d, Uf + \ {d, + div., (2.7) 

where div. denotes an angular divergence, (sin 6)~^ 96i(sin 9D^)+d^D'^ . Instead of the asymp- 
totic wave forms /, /+, /x (which have the dimension of length), one can work with the 
corresponding "news functions" 

N{v, e, ^) = f{v, e, v>) , N+ = h, N,= /x , (2.8) 



^These coordinates must be restricted by the condition that the incoming coordinate cones v = 
t + r = const, be asymptotically tangent to exact cones of g^y{x'^). 

^The three functions /, /+, /x of v, 6, cp are equivalent to six functions in M^, i.e. six functions of 
r, 9, ip with r > 0, because the advanced time v ranges over the full line (— oo, +oo). 
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which are dimensionless, and whose squares give directly the incoming energy flux appearing 
on the R.H.S. of Eq. (^). 

Before specializing this generic ingoing state, let us note |T^, the existence of two impor- 
tant global symmetries of the classical field equations, (|2.3|) and ( p.4|) . They are invariant 
both under global scale transformations, and under a constant shift of : ds'"^ = a^ds"^, 
(f)' = (f) + b. [They are also invariant under local coordinate transformations.] In terms of 
asymptotic data (whose definition requires a specific "flat" normalization of the coordinates 
at past infinity), the global symmetry transformations read: 

f{v, 9, if) ^ f\v', e,ip)=a f{v'/a, 9, if) , (2.9) 
U{v, 9, ^) ^ r^{v\ 9,^) = a Uiv'/a, 9, ^) , (2.10) 
{v, 9, ^) {v', 9,^) = a /x (v'/a, 9, ^) (2.11) 

with r — > r' = ar, v —>■ v' = av, and 

(po ^ (p'o = (Po + b . (2.12) 

Note that the three dimensionless news functions (|2.8|) are numerically invariant under the 
scaling transformations: N'{v',9,{p) = N{v,9,ip), etc. with v' = av. 

As the full time evolution of string vacua depends only on the null data, it is a priori 
clear that the amplitude of the dimensionless news functions must play a crucial role in 
distinguishing "weak" incoming fields -that finally disperse again as weak outgoing waves- 
from "strong" incoming fields that undergo a gravitational instability. To study this difficult 
problem in more detail we turn in the next section to a simpler model with fewer degrees of 
freedom. 

Before doing so, let us comment on the "classicality" of an "in state" defined by some news 
functions (|2.8|) . The classical fields ( p.5|) , (|2.6| ) deviate less and less from a trivial background 
{(pOyVfj-i^) as r — > oo (with fixed v,9 and (f). One might then worry that, sufficiently near 
X~ , quantum fluctuations 0"^, hj^^, might start to dominate over the classical deviations 
0'^ = 0(x) — 00) h'^^u = df^i' ~ Vf^u- To ease this worry let us compute the mean number of 
quanta in the incoming waves ( |2.5| ), ( p.6| ). Consider for simplicity the case of an incoming 
spherically symmetric dilatonic wave 0?^ = f{v)/r + o(l/r). Such an incoming wave has 
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evidently the same mean number of quanta as the corresponding outgoing]^ wave (pouti with 
asymptotic behaviour 0^^^ = —f{u)/r + o(l/r) , on X+ (r oo, with u = t ~ r,6 and ip 



fixed). As is well known ||21]], a classical outgoing wave can be viewed, quantum mechanically, 
as a coherent state in the incoming Fock space defined by the oscillator decomposition of the 
quantum field 0. Instead of parametrizing 0^^^. by the waveform f{u), one can parametrize 
it by an effective source J{x^) such that 0out(^) — I d^^' Pret{x — x') J(x'), where P^^t is the 
(action-normalized) retarded propagator of the field 0. [We use a normalization such that 
S = — I (f) P^^ (f) + J0.] The mean number of quanta in the coherent state associated to a 



classical source J{x) has been computed in for the electromagnetic case, and in p2[ for 



generic massless fields. It reads (with px = px — ut) 

n=lJ ^Jip')Ji-p)nJip), J{p) ^ J d'xe-'^\J{x), (2.13) 



where 71 denotes the residue of the propagator of the corresponding field, defined such that 
the field equation P~^ = J reads ^0 = —TZ J. [In the normalization of the present paper 
7^ = 167rG.] Inserting J{x) = -{AG)-^ f{t) S^^\x) into Eq. (|27[3|) yields 

^=2GEi 2^"l^(-)l =2Ga/o 2^;;I'^(-)I ■ 

where N{uj) = J dte^^^ N{t) is the one-dimensional Fourier transform of the news function 
N{t) = f'{t). In order of magnitude, if we consider a pulse-like waveform f{t), with charac- 
teristic duration At ~ 4, and characteristic dimensionless news amplitude Nc ~ fc/^c, the 
corresponding mean number of quanta is 

---^ = ^- (2-15) 

In the present paper, we shall consider incoming news functions with amplitude ^ 1 and 
scale of variation ^s- Moreover, in order to ensure a sufficient amount of infiation, the 
initial value of the string coupling g"^ = e'^° has to be <^ 1. Therefore, Eq. ( p. 151) shows 



^In keeping with usual physical intuition, it is more convenient to work with an auxiliary problem 
where a classical source generates an outgoing wave which superposes onto the quantum fluctuations 
of an incoming vacuum state. 
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that such incoming states can be viewed as highly classical coherent states, made of an 
enormous number of elementary 0-quanta: n 1. This quantitative fact clearly shows 
that there is no need to worry about the effect of quantum fluctuations, even near X~ where 
'Pin ~ /('^)/^ ~^ 0- the case of the negative- curvature Friedmann-dilaton solution this 
fact has been recently confirmed by an explicit calculation of scalar and tensor quantum 
fluctuations |^ . 



III. EQUATIONS FOR A SPHERICALLY SYMMETRIC PRE-BIG BANG 

The simplest nontrivial model which can exhibit the transition between weak, stable data, 
and gravitationally unstable ones is the spherically-symmetric Einstein-dilaton system. The 
null wave data for this system comprise only an angle-independent asymptotic dilatonic wave 
form /(f), corresponding to the dimensionless news function 

N{v) = d,f{v). (3.1) 

The spherically symmetric Einstein-dilaton system has been thoroughly studied by many 



authors with the strongest analytic results appearing in papers by Christodoulou [Q- 
A convenient system of coordinates is the double null system, (m, f), such that 

ct) = (t){u,v), (3.2) 



ds^ = —Q'^{u, v) du dv + r^(M, v) duj'^ , (3.3) 

where duj"^ = dO"^ + s\.v?6d({P'. The field equations are conveniently re-expressed in terms 
of the three functions (f){u,v), r{u,v) and m{u,v), where the local mass function m{u,v) is 
defined by | 

One gets the following set of evolution equations for m, r and 



In order to avoid confusion we indicate, in standard notation, the variable kept fixed under 
differentiation. 
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d^r 2m /dr\ /dr\ {3 7) 

dudv r — 2m \du J _^ \dv J ^ ' 

r-^+(—] (^] +(—] (^] -0 (3 8) 

dudv \du J _^\dv J \dv J ^\du J 

This double-null evolution system is form-invariant under independent local reparametriza- 
tions of u and v. u' = U{u), v' = V{v). To freeze this gauge freedom it is convenient to 
require that m = f on the central worldline of symmetry r = (as long as it is a regular part 
of space-time), and that v and r asymptotically behave like ingoing Bondi coordinates on 

(see below). 

The quantity 

_ 2m{u,v) 

H{u,v) = (3.9) 



plays a crucial role in the problem. Following we shall sometimes call it the "mass 
ratio". If /i stays everywhere below 1, the field configuration will not collapse but will finally 
disperse itself as outgoing waves at infinity. By contrast, if the mass ratio /i can reach 
anywhere the value 1, this signals the formation of an apparent horizon A. The location of 
the apparent horizon is indeed defined by the equation 

A: fi{u,v) = l. (3.10) 

The above statements are substantiated by some rigorous inequalities stating that: 

/ 3t \ ( Stti \ 

- (l-/x)>0, — (l-/i)<0. (3.12) 



dv I \ du 



V 



Thus, in weak-field regions (/i < 1), (9t,r)„ > 0, while, as /i > 1, (9t,r)„ < 0, meaning that the 
outgoing radial null rays ("photons") emitted by the sphere r = const, become convergent, 
instead of having their usual behaviour. 

It has been shown long ago that the presence of trapped surfaces implies the existence of 
some (possibly weak) type of geometric singularity [^. In the case of the spherically 
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weak-field expansion 
cosmological regime 



FIG. 3. Schematic representation of the space-time generated by the collapse of a spherically 
symmetric pulse of dilatonic waves. An incoming pulse of scalar news N (v) , which grows by 0{1) on 
an advanced-time scale li, collapses to a space-like singularity B after having formed an apparent 
horizon A hidden behind the event horizon 7i. In the vicinity of A and 7i, there is an abrupt 
transition between a weak-field region (dark shading) where the perturbation series J[)~( P^ 
holds, and a strong-field one (light shading) where the cosmological-like expansion ( |5.34 )~( l536 } 
holds. 



symmetric Einstein-dilaton system, it has been possible to prove ||2^ that the presence of 
trapped surfaces (i.e. of an apparent horizon A ( |3.10|) , the boundary of the trapped region) 
imphes the existence of a future space-hke singular boundary B of space-time where the 
curvature blows up. Both A and B are "invisible" from future null infinity (X+), being 
hidden behind an event horizon 7i (a null hyper surf ace). See Fig. |^. 

One of our main purposes in this work is to give the conditions that the incoming dilatonic 
news function N{v) must satisfy in order to create an apparent horizon, and thereby to lead 
to some localized gravitational collapse. Before addressing this problem let us complete the 
description of the toy model by re-expressing it in terms of the ingoing Bondi coordinates 
(w, r). The metric can now be written in the form 



(3.13) 



and the field equations ( |3.5| )- (|37^ ) become 
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dr 



4 \dr 



' dm\ / 2m\ ( d(j) 



dr 



r J 4 \dr 



(3.14) 
(3.15) 



2 = e-'^ — 

\dv 2 



2m\ (d<i)\ (d<i)^ 
r J \ dr I \ dv , 



(3.16) 



9V 1 fd(l)' 



+ 



dvdr r \dv I 
Eqs. ( CT ), ( PJ[50 imply: 



+ 6^1 



r J \ or'^ / r V r J \ or 



2m\ 



(3.17) 



a, (e'^("''') (r - 2m)) = e^^'''") . 



(3.18) 



In this coordinate system one can solve for the two functions /3(f , r) and /x(f , r) (or m(f , r)) 
by quadratures in terms of the third unknown function 0(f,r). Indeed, imposing that the 
coordinate system be asymptotically flat at I~ (i.e. that /3(f , r) — > as r — > cxd), one first 
finds from Eq. ( |3.14| ) 

n+OO 



where 



(3{v,r) = - / —F{v,r'), 

Jr r 



(3.19) 



(3.20) 



Eq. ( |3.15D can then be integrated to give 

'2m^{v) - /+°° dr'F{v, r') e^^"^"'^ 



H{v,r) 



where the "integration constant' 



r €■ 



P(v,r) 



mooiv) = lim miv.r) 



(3.21) 



(3.22) 



denotes the incoming Bondi mass. By the Bondi energy-flux formula ( |2.7| ), which is also the 
limiting form of (p.l6| ) for r oo, moo(f) is given in terms of the news function by 



(3.23) 
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where we have inserted the asymptotic behaviour of on X 



(^,r) = 0o + ^ + of-) , (3.24) 



and assumed that N{v) decays sufficiently fast as f — oo (i.e. near past time-hke infin- 
ity). From Eqs. ( p. 191 ) and ( p.21[ ) we find the leading asymptotic behaviour of the metric 
coefficients to be 

r r \r J 

^I3(v,r) ^i_^Qf}_y (3 26) 



IV. DATA STRENGTH CRITERION FOR GRAVITATIONAL INSTABILITY 

The purpose of this section is to outline the condition that the initial data, i.e. the 
wave form /(f) or the news function N{v) = f'{v), must satisfy in order to undergo, or 
not undergo, gravitational collapse. A lot of mathematical work has been done on this issue 
2^- p6| . In particular, Refs. |2^ and gave two different "no collapse" criteria ensuring 



that the field configuration never collapses and finally disperses out as weak outgoing waves 



if some functional of the data is small enough. On the other hand, Ref. gave a sufficient 
"collapse" criterion ensuring that the field configuration will collapse, if some functional of 
the data is large enough, thereby giving birth to a curvature singularity hidden behind a 
horizon. The problem with these nice and important mathematical results is threefold: (i) 
these criteria are sufficient but not necessary, so that they cannot answer our problem of 
finding (if possible) a sharp criterion distinguishing weak-field from strong-field data; (ii) 



the various measures of the "strength" of the data given in Refs. , |^ and ||2^ are quite 
unrelated to each other and do not point clearly toward any sharp "strength criterion" ; and 
(iii) these criteria are not expressed in terms of the asymptotic null data N{v). 



Our aim here is not to compete with Refs. [£4|, pSf, on the grounds of mathematical 



rigour, but to complement these results by a non rigorous study which leads to the iterative 
computation of a sharp strength criterion, i.e. a functional S[f], such that data satisfying 
S < 1 finally disperse out at infinity, while data satisfying S > 1 are gravitationally unstable 
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and partly collapse to form a singularity. We will then compare our strength functional S 
with the rigorous, but less explicit results of Note that the emphasis here is not on the 
question of whether the created singularity is hidden beyond a global event horizon when 
seen from outside, but rather on the quantitative criterion ensuring that past trivial data 
give birth to a space-like singularity which we shall later interpret as a budding pre-Big Bang 
cosmological universe. 



A. Perturbative analysis in the weak-field region 

In view of the discussion of the previous section, the transition between "weak" and 
"strong" fields is sharply signaled by the appearance of an apparent horizon, i.e. of space- 
time points where fi = 2m/ r reaches the "critical" value 1. Therefore, starting from the 
weak-field region near X~, we can define the strength functional by computing fi{u,v) as a 
functional of the null incoming data and by studying if and when it can exceed the critical 
value. Specifically, we can set up a perturbation analysis of the Einstein-dilaton system in 
some weak- field domain connected to I~ , and see for what data it suggests that /i will exceed 
somewhere the value 1. The expected domain of validity of such a weak-field expansion is 
sketched in Fig. ^ 

The {u,v) system, Eqs. ( |3.2| ) and (|3.3|), lends itself easily to a perturbative treatment in 
the strength of 0. We think here of introducing a formal parameter, say A, in the asymptotic 
data by f{v) — > \f{v) and of constructing the full solution (f){u,v), r{u,v), m{u,v), as a 
(formal) power series in A (or by some better iteration method): 

0(m, v) = (f){u, i;) - 00 = A 01 + A303 + A^05 + . . . , (4.1) 
r{u, v) = ^{v ~ u) + r2 + n + . . . , (4.2) 
m{u, v) = X^ m2 + A^ m4 + . . . . (4.3) 

Indeed, knowing to order A^^^^ and r to order A^^~^ we can use Eq. (|3.5| ) or (|3.6| ) to 
compute by quadrature m to order A^^. We can then rewrite Eqs. ( ^^ , ( ^TsD as 

duv{r 4>) = u r (p , duvr = ujr, (4.4) 

where 
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uj = ^ — , 4.5 

1 — /i 



and therefore the right-hand sides of Eqs. (|4.4| ) are known to (relative) order A^^. Solving 
Eq. ( [4.4D , then yields and r to the next order in A (A^^^^ and A^^, respectively). To be 
more explicit, Eqs. ( |4.4| ) are of the form 

duv ip{u, v) = a{u, v) , (4.6) 

where the "source" a = ujip is known at each order in A from the lower-order expression 
of /i and ip, and where the "field" iIj{u,v) is restricted to vanish on the central worldline 
of symmetry v = u (where r = 0, and 4> is regular in the weak-field domain). The generic 
Eq. ([4. 61) can then be solved by introducing the retarded Green function 



G{u, v; u', v') = e{u - u') e{v - v') - e{v - u') e{u - v') , (4.7) 



duv G{u, v; u', v') = S{u — u) S{v — v') . (4.8) 

This is the unique Green function in the domain v > u, v' > u' which vanishes at v = u, 
and whose support in the source point {u',v') lies in the past of the field point {u,v). The 
general solution of (|4.6| ) can then be written as 

■ijj{u,v) = ipiniuju) + / / du dv' G{u,v;u ,v') a{u ,v') , (4.9) 

J Jv'>u' 

where ipin{u, v) is the free "incoming" field, satisfying d^v ipm = and vanishing when v = u, 
i.e. 

tpiniu, v) = g;^{v) - gi^{u) . (4.10) 

The incoming waveform 5'in(f ) is uniquely fixed by the incoming data on X^. For instance, 
for r(u, v) asymptotic flatness gives rin(M, v) = ]^{y — u) (we recall that in fiat space v = t + r 
and u = t — r), while for (j){u, v) = (f){u, v) — (pQ the asymptotic expansion ( p.24| ) (in which 
f{v) was supposed to vanish as f ^ — oo) yields (t>i^{u,v) = 2{f{v) — f{u))/{y — u). 

In principle, this perturbative algorithm allows one to compute the mass ratio fi{u,v) = 
2m/ r to any order in A: 

/i(M,f) = AV2 + AV4 + ... • (4.11) 
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Evidently, the convergence of this series becomes very doubtful when /i can reach values of 
order unity, which is precisely when an apparent horizon is formed. One would need some re- 
summation technique to better locate the apparent horizon condition fi{v,u) = 1. However, 
we find interesting to have, in principle, a way of explicitly computing successive approxima- 
tions to the possible location of an apparent horizon, starting only from the incoming wave 
data /(f). 



B. Strength criterion at quadratic order 



Let us compute explicitly the lowest-order approximation to ( [4.11| ) (we henceforth set 



A = 1 for simplicity). It is obtained by inserting the zeroth-order result ro(u, v) 



and the first-order one for 



fiv)-fiu) 2[f{v)-f{u)] 



To 



V — u 



into Eq. (^ or (gj). This yields (with f{v) = df/dv 

\ 2 



dm2 

dv 



1 

° ^ 2 



rn I - 



f'iv) - 



fiv) - fiu) 



V — u 



[V — u) 



(4.12) 



(4.13) 



' dlTLQ 



'1 901 



fiv)-fiu) 



\ du J ^ ' ^ \2 du J 4[ ' V — u 

The compatibility of the two equations is easily checked, e.g. by using 

9.(r2(9,0i)') =ro 9,01 9,01 = -9,(r2(9,0i)2), 



(4.14) 



(4.15) 



which expresses the "conservation" of the 0-energy tensor in the {u, v) plane. Noting that 
m{u, v) must vanish (by regularity of r = in the weak-field domain) at the center v = u 
one gets, by quadrature, the following explicit result for 



_ 1-2 

m2{u, v) = ■m{u, i;) - - ro 01 , 



(4.16) 



where 



mm, V) 



(4.17) 
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Note that Eq. ( [4.16| ) decomposes m2{u, v) in a "conserved" piece m{u, v) = g{y) ~g{u) (with 
duv"^ = 0) and a ^-dependent one. This is an integrated form of the local "conservation" law 
( [4.15| ) of the 0-energy tensor. Using ( ^.16| ) one gets a very simple result for the quadratic 
approximation to the mass ratio /i = 2m/r, namely 

2 



V — u 



dx f (x 



V — u 



(4.18) 



Let us introduce a simple notation for the average of any function g{x) over the interval 

1 



{9)[u,v\ = {9{x))x&[u,v] 



V — u 



dx g{x) . 



(4.19) 



Then the mass ratio, at quadratic order, can be simply written in terms of the scalar news 
function A^(t>) = /'(f): 



/X2(M,t^) = (iV')K.]-((iV)[n,.])', 

= ((iV(x)- (iV)[„,,])2),e[„,„], 
= Var iN{x)),x(.[u,v] , 



(4.20) 
(4.21) 
(4.22) 



where Var {g)[u,v] denotes the "variance" of the function g{x) over the interval [u,v], i.e. the 
average squared deviation from the mean. 

As indicated above, having obtained fi to second order, say /i2 = V{u,v) = Var (A^) 
one can then proceed to compute r and to higher orders. Namely, from Eq. ([4.9| ) 

'V{u,v) 



2r{u, v) = v — u — G* 



V — u 



(4.23) 



[rct)){u,v) = f{v)-f{u)-G* 



f{v) - f{u) 



+ C>(A^ 



(4.24) 



{v — uY 

where the star denotes a convolution. Then, one can get from these results m4, and thereby 
/i4, by simple quadratures. In the following we shall only use the quadratic approximation to 
yU, though we are aware that when the initial data are strong enough to become gravitationally 
unstable, higher order contributions to /i are probably comparable to /i2- 

Finally, at quadratic order, we can define the strength of some initial data simply as the 
supremum of the variance of the news over arbitrary intervals: 



19 



^2 = sup/i2 = sup VaT{N{x))^^[u,v] ■ (4.25) 



u,v 



At this order ^2 < 1 means that fi2{u,v) stays always below one, i.e. that no trapped region 
is created, and that the field is expected to disperse, while 1S2 > 1 signals the formation, 
somewhere, of trapped spheres, and therefore (by the results of [^) the formation of a 
singularity. To be more exact we should actually define, in analogy with ([4.25|) , a quantity 
S2p = sup X]fc=i /^2fc, at any finite order in the weak field expansion, and state our collapse 
criterion as the inequality: 

S2p>C2p, C2/^°°l. (4.26) 

In practice, even if we do not expect it to be quantitatively exact, we will use the criterion 
( [4.26| ) with p = I, taking C2 = 1, hoping that it will be quahtatively correct in capturing the 



features of the news function which are generically important for producing a gravitational 
collapse. Let us note that the functional S2 is (like A^) dimensionless (and therefore scale 
invariant) and that it is nonlocal. We note also that it is invariant under a constant shift of 
A^, N{v) N{v) +a, which corresponds to adding a linear drift in / : f{v) —>■ f{v) +av + b. 
Such a shift is (formally) equivalent, in view of Eq. (|4.12| ), to a constant shift of 0o 0o + 2a. 



The non locality of ^2 is physically interesting because it indicates that it is not the 
instantaneous level of the energy fiux N'^{v) which really matters, but rather the possibility of 
having a fiux which varies by ~ 100% over some interval of advanced time. This non locality 
defines also some characteristic scales associated with the collapse (when ^2 > 1). Indeed, 
if, by causality, we consider increasing values of v, and define S2{v) = sup Var(A^)[„ t,], the 

u,u<v 

first value of v for which S2{v) exceeds 1 defines an advanced characteristic time when the 
collapse occurs, and the corresponding maximizing interval [u, v] defines a characteristic time 
scale of (in)homogeneity. [Note that fi2{u,v), Eq. ( |4.22| ), vanishes when u = v, and vanishes 
also generally when u —00 if N{u) tends to a limit at —00.] 

Let us also mention that our strength functional (|4.25| ) is superficially similar to the 



one recently introduced by Christodoulou p6| to characterize sufficiently weak (i.e. non 



collapsing) data. In the lowest order approximation his criterion measures the strength of 
the data by V = TV{N{x)), where TV denotes the total variation (i.e. essentially V = 
dv\N' {v)\) . Like our variance, this is a measure of the variation of N{v) with, however, 
a crucial difference. If V were a good criterion for measuring the strength of possibly strong 
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data, we would conclude that a news function which oscillates with a very small amplitude 
for a very long time will be gravitationally unstable, while our criterion indicates that it will 
not, which seems physically more plausible. We hope that our strength functional ( ^.25| ) will 
suggest new gravitational stability theorems to mathematicians. 



C. Comparison with a collapse criterion of Christodoulou 



To check the reasonableness of our quadratic strength criterion ( [4.25|) , ([4.26|) we have 
compared it with the rigorous, but only sufficient, collapse criterion of Christodoulou 
and we have applied it to two simple exact solutions. Ref. |^ gives the following sufficient 
criterion on the strength of characteristic data considered at some finite retarded time u 



2Am 
Ar 



> 



ri f ri 
V, [2A-r 



+ ^-1 
r2 



(4.27) 



where vi < r2, r2 < 3ri/2, are two spheres, Ar = r2 — ri is the width of the "annular" region 
between the two spheres, and Am = m2 — nii = m{u, r2) — m{u, ri) is the mass "contained" 
between the two spheres, i.e. more precisely the energy flux through the outgoing null cone 
u = const., between ri and r2. Explicitly, from Eq. (|3.15|) in (m, r) coordinates, we have 

„2 - ^ . /o,n2 



2Am 



r2 / 2m 

^^T ^ / \ 

r-i 4 V r J \or 



(4.28) 



We can approximately express this criterion in terms of the incoming null data N{v) = f'{v) 
if we assume that the outgoing cone u = const, is in the weak-fleld domain, so that we can 
replace on the R.H.S. of Eq. ( CT) by = r-^[fiu + 2r) - f{u)]. Then the L.H.S. of 
Eq. ( [4.27|) becomes (at quadratic order) 

2Am M{u,V2) - M{u,vi) 



Ar 



V2 - Vi 



(4.29) 



where 



M{u,v) = [v - u)VaiiN)[u,v] 



dx 



X] 



(4.30) 



Therefore, in this approximation, the criterion of Ref. |^ becomes (with vq = u) 

M{vo,V2) - M{vo,Vi) 



3vo,Vi,V2 ■■ 



V2 - Vi 



> C{vo,Vi,V2) 



(4.3i; 
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with the constraints 2v2 — 3?;i < — fo, Vi < V2, and the definition 

C{vo,Vi,V2) = log 777 7 + 1- (4-32) 

V2-V0 \2{V2-Vi)J V2-V0 

If vq is at the boundary of its allowed domain, i.e. if vq = 3vi — 2v2, this criterion is fully 
compatible with ours since C(fo,fi,f2) = 3 implying through Eq. ( [4.30|) 



Var(iV)[,„,,,] > 1 + ^ Var(iV)[,„,,,] > 1 . (4.33) 



In the opposite case (f large and negative) we think that the criterion (|4.31| ) is also com- 



patible, for generic news functions, with the general form suggested above, i.e. 

3v[,v'2 : Var(iV)[,,,,,] = (iV^);,,^,,] - {N)f,,^.,,^ > C2 , (4.34) 

with some positive constant C2 of order unity. We first note that one should probably impose 
the physically reasonable condition that N{v) decays [] sufficiently fast as f ^ —00 (faster 
than |f |~^^^) to ensure that the integrated incoming energy fiux J^'-^ dv N'^{v) is finite. [This 
constraint freezes the freedom to shift N{v) by a constant.] 

When — fo is large, the question remains, however, to know how fast the function N{v) 
decays when vi—v (for vq < v < vi) becomes much larger than V2—V1. Let us first consider 
the physically generic case where N{v) decays in a reasonably fast manner, say faster than 
a power |f I"'' with k, of order unity. Then choosing fo iii the criterion ( [4. 311 ) just large 



enough to allow one to neglect A^(fo) with respect to N{v), say (f 1 — fo) ~ (f2 — vi) x 10^^'^, 
corresponding to A^(fo) = 0.1 A^(fi), the L.H.S. of Eq. ( ^4.31|) becomes approximately 

^-(iV^) (4.35) 

while the function C on the R.H.S., which grows only logarithmically with fo, becomes 

CK.„„.)^l„g(^^;LZ^)+5. (4.36) 



^ Actually, when comparing our criterion (based on asymptotic null data) to that of Ref. 
(based on characteristic data taken at a finite retarded time u = vq) we can consider, without loss 
of generality, that the news function N{v) vanishes identically for v < vq. 
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Therefore, a function of the type N^{v) ~ c [{v2 — v)/{v2 — "Wi)] " (when t> < t>i < V2) with 
> 5 + (log 10) / K, will fullfil the criterion (|4.31|) . To see whether this is compatible with our 
criterion ( [4.34|) we have studied (analytically and numerically) the variance, over arbitrary 
intervals [v[, v'2], of the function N^{v). We found that the inequality ( [4.34| ) is satisfied for a 
constant which is of order unity if k stays of order unity. [For instance, in the extreme 
case K = 1/2 corresponding to a logarithmically divergent incoming energy flux, we find 
02^^"^^ ~ 0.2]. We note, however, that if one considers extremely slow decays of N{v), i.e. 
very small exponents k, (completed by a faster decay, or an exact vanishing, before some 
large negative cutoff Vc ), the constant 6*2'^^ needed on the R.H.S. of the variance criterion 
i ^M ) tends to zero. This signals a limitation of applicability of our simple variance criterion 
based only on the quadratic-order approximation /i = /i2 + /i4 + ■ ■ ■ — /i2- Indeed, one can 
check that, in such an extreme situation /i2 is abnormally cancelled, while /i4 will be of order 
unity. However, we believe that for generic, non extremely slowly varying news functions, 
the simple "rule of thumb" ( |4.34|) is a reasonable approximation to the (unknown) exact 
collapse criterion. To further check the reasonableness of our criterion ([4.25|) we turn our 
attention to some exact solutions. 



D. Exact solutions 



A first exact, dynamical Einstein-dilaton solution is defined by a news function consisting 
of the simple step function 

N{v)=p9{v). (4.37) 

The corresponding solution has been independently derived by many authors p9|. Here, p 
is a real parameter and the value |p| = 1 defines the threshold for gravitational instability: 
no singularity occurs for \p\ < 1, while \p\ > 1 causes the birth of singularity. The metric for 
this solution takes the form (note that Q = 1) 

ds^ = -dudv + duo^ , (4.38) 

with the solutions for 0, r and m given as follows: 

= 0, r = ]-{v-u), m = 0, for t; < , (4.39) 
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while, for f > 0, 

(j){u,v) = -lo 
r^(u, v) 
m{u, v) 



1 
4 



{1 — p) V — u 
{1 + p)v — u 

[1 — p^) ~ 2vu + 



[(1 — p)v — u\ [(1 + p)v — u] 



p uv 



p^ uv 



[(1 — p^) v"^ — 2vu + u"^] 
Let us note that the perturbation-theory value of /i for > 0, is 



p^ uv 
(v — u)' 

-.2 



-u] 



(4.40) 
(4.41) 
(4.42) 
(4.43) 

(4.44) 



The maximum value of /i2, reached for u = — f, is p /A. This shows that, as expected, 
the quadratic order criterion ([4.25|) is only valid as an order of magnitude, but is quantita- 
tively modified by higher-order corrections. This example, and a related general theorem of 
Christodoulou [^, suggest that, if the exact criterion were of the type ^2 = C2, the constant 
C2 should equal ^. (Note that this value is also compatible with the constant c'^^'^'^ ~ 0.2 
appropriate to l^;!^^/^ decay.) 

A second exact solution |^|, |T^ is a negative-curvature Friedmann-like homogeneous 
universe. It is defined by the following null data on (for f < 0) 



/(.) ^ 



V 



(4.45) 



Here we view this solution as defined by incoming wave data in a fiat Minkowski background. 
Actually the data are regular only in some advanced cone f = T + i? < 0, and blow up when 
f = 0. Here, T and R = VA^ + l^ + Z^ are usual Minkowski- like coordinates in the 
asymptotic past. In terms of such coordinates the exact solution reads (when — i?^ < 1) 



1 - 



1 



(T2 - i?2)2 



l-dT^ + dX^ + dY^ + dZ^ 



-VS l0£ 



T2 - i?2 



1 



J^2 



(4.46) 



(4.47) 



Note that, for simplicity, we have set here 0o = 0, and we have also set the length scale 
appearing in f{v) = —^/Sil/v to Iq = 1. We shall come back later to the cosmological 
significance of this solution. Let us only note here that, in terms of the null coordinates 
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u = T-R, v = T + R, (4.48) 

the exact mass ratio reads 

Ku,v)= \ ^ , 4.49 
while, starting from ( [4.45| ), one obtains 



^i2{u,v) = Var(iV)[„,,] = ^ / . (4.50) 

One finds that a strength criterion of the form sup fi2{u,v) = C2 is first satisfied (as one 
increases v from ~oo) when v = v^:2, with u = m^,2, where 

'y*2 = -4 (^^y^ j ' = 3t;=^2 • (4.51) 

For C2 of order 1, this is in quahtative agreement with the exact result that the apparent 
horizon is located at 

^A=\, r^ = -^^{l-vy/\ (4.52) 

so that the apparent horizon and the singularity are first "seen" at exact = —^o- 



V. TRANSITION FROM THE WEAK-FIELD TO THE COSMOLOGICAL 

REGIME 

The following general picture emerges from the previous sections: Let us consider as 
"in state" a generic classical string vacuum, which can be described as a superposition 
of incoming wave packets of gravitational and dilatonic fields. This "in state" can be 
nicely parametrized by three asymptotic ingoing, dimensionless news functions N{v,6,ip), 
N^{v, 6, (f), Ny.{v, 6, if). When all the news stay always significantly below 1, this "in state" 
will evolve into a similar trivial "out state" made of outgoing wave packets. On the other 
hand, when the news functions reach values of order 1, and more precisely when some global 
measure of the variation ^ of the news functions, similar to the variance ( [4.25[ ), exceeds some 



*The argument that the collapse criterion should be (at least) invariant under constant shifts of 
all the news functions (corresponding to classically trivial shifts of the background fields g^^^, (po) 
indicates that only the global changes of the news matter. 
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critical value of order unity, the "in state" will become gravitationally unstable during its 
evolution and will give birth to one or several black holes, i.e. one or several singularities 
hidden behind outgoing null surfaces (event horizons). Seen from the outside of these black 
holes, the "out" string vacuum will finally look, like the "in" one, as a superposition of 
outgoing waves. However, the story is very different if we look inside these black holes and 
shift back to the physically more appropriate string conformal frame. First, we note that 
the structure of black hole singularities in Einstein's theory (with matter satisfying p < e) is 
a matter of debate. The work of Belinsky, Khalatnikov and Lifshitz has suggested the 
generic appearance of an oscillating space-like singularity. However, the consistency of this 
picture is unclear as the infinitely many oscillations keep space being as curved (and "turbu- 
lent" [^], [Q) as time. Happily, the basic gravitational sector of string theory is generically 
consistent with a much simpler picture. Indeed, it has been proven long ago by Belinsky 



and Khalatnikov [34| that adding a massless scalar field (which can be thought of as adding 
matter with p = e as equation of state) drastically alters the BKL solution by ultimately 
quenching the oscillatory behaviour to end up with a much simpler, monotonic approach 
to a space-like singularity. When described in the string frame, the Einstein-frame collapse 
towards a space-like singularity will represent (if grows toward the singularity) a super- 
inflationary expansion of space. The picture is therefore that inside each black hole, the 
regions near the singularity where grows will blister off the initial trivial vacuum as many 
separate pre-Big Bangs. These inflating patches are surrounded by non-inflating, or deflating 
(decreases 0) patches, and therefore globally look approximately closed Friedmann-Lemaitre 
hot universes. This picture is sketched in Fig. |I] and Fig. |^. We expect such quasi-closed 
universes to recollapse in a finite, though very long, time (which is consistent with the fact 
that, seen from the outside, the black holes therein contained must evaporate in a finite 
time). To firm up this picture let us study in detail the appearance of the singularity in the 
simple toy model of Section |T|. 



A. Negative-curvature Friedmann-dilaton solution 

We have already seen that in the toy model there is an infinitely extended incoming 
region where the fields are weak and can be described as perturbations upon the incoming 
background values ((7°^ = rj^u, 4>o)- One expects that the perturbation algorithm described in 
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section P^V A| becomes unreliable when corrections become of order unity. In particular, one 
generically expects that the apparent horizon /i = ^ = 1 will roughly divide space-time in 
two domains: the perturbative incoming weak-field domain where fi{u, v) <^ 1, and a strong- 
field domain where fi{u,v) ^ 1. [Actually, we shall see below that the precise boundary of 
the domain of validity of perturbation theory may also depend on other quantities than just 
/i.] This separation in two domains is represented in Fig. ^ 

For guidance let us study in detail this separation in two domains in the case of the 



k = —1 Friedmann-dilaton solution [Q, |T^. Let us first write explicitly the solution Eqs. 
( [1.46| ) and ([4.47|) in {u, v) coordinates 

'uv — V 



-Vt^{u,v) dudv + r'^ du^ , = — -\/3 log 



uv + 1 



r^(M, v) 



1 



1 



h) 



uv{v — uY 



1 - 



2m 



4 V u 
{u^v — 1) {uv^ — 1) 



The perturbative algorithm of section p.V A gives 



2^3 



uv 



3 

i-i 1 



r{u, t>) ~ - (t> 



u 



2 f 2 



(5.1) 

(5.2) 
(5.3) 

(5.4) 
(5.5) 



which agree with the expansions of Eqs. (5^) and (5^). One sees that perturbation theory is 
numerically valid up to, say, uv < -\/2, at which point there is an abrupt transition towards 
the cosmo logical singularity located at uv = 1. Though the transition surface uv ~ \/2 
globally differs from the apparent horizon fi{u, f ) = 1, we note that our criterion points out 
to a specific event (m*2, ^^2) on the apparent horizon A (the point where A is first "seen" 
from infinity) which lies, roughly, at the intersection of A and of the transition surface. This 
confirms that our criterion is able, at least in order of magnitude, to correctly pinpoint when 
and where one should shift from the perturbative regime to a different, cosmological-type 
description. In the case of the solution Eqs. ( |5.1| )-(|^), one sees better the cosmological 
nature of the strong-field domain by introducing the following coordinates 



u 



T-R. 



T + R. 



(5.6) 



satisfying the useful relations 
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— T = e cosh ^ , 



R 



sinh ^ , 



5« 



{dr, - dt 



e '^^ = uv . 



,2? 



e-' ' ^ (9^ + di:) , 



In terms of the coordinates (?7,0 ^^e solution reads ||3^, [|18 



ds'^ = a^{ri) -dTf + + sinh^ ^ duo' 
a^{r]) = 4 cosh sinh(— 77) , e"*^ = 



sinh ?7 ' 
cosh?7 



V3 



(5.7) 
(5.8) 
(5.9) 



(5.10) 
(5.11) 



We note that the solution is regular in the domain —00 <r7<0, 0<.^< +00 (which 
corresponds to the past of the hyperboloid — i?^ = 1 in Minkowski-like coordinates) and 
that a space-like singularity is reached at = (i.e. uv = T"^ — = 1). Let us also note 
the expressions 



= ol^iv) sinh^^ = 4 cosh?] sinh(— r^) sinh^^ 



(5.12) 



1 - 



2m 



sinh^ ^ coth^ ^ — coth^ 2// 



(5.13) 



telling us that the apparent horizon is located at 21] = Near the singularity, rj 0^ , the 
dilaton blows up logarithmically while the metric coefficients have the following power-law 
behaviours: 



— log(— 77) , Q^{C,,ri) = a^{r]) —Arj , r^(^, r/) ~ — 4// sinh^ ^ 



(5.14) 



Before leaving this example we wish to emphasize some features of it which appear to follow 
from its homogeneity and could be misleading for the general case. The singular boundary 
terminates, in this case, on X~, instead of on X"*", as one generically expects. The apparent 
reason for this is the singularity in the flux at i; = which creates a future boundary on 
X~. We believe that, in this case, a better description of physics is obtained by going to 
new non-Minkowskian coordinates of Milne's type (see |12|) which automatically incorporate 
the singularity on X~ . However, if one restricts oneself to more regular initial data, having 
a finite integrated energy fiux (generalized pulse-like data), the singular boundary should 
never come back to X^. [At most, it could end at space- like infinity for non integrable 
total energy fiux]. 
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B. Kasner-like behaviour of Einstein-dilaton singularities 



Coming back to the generic case of an arbitrary Einstein-dilaton cosmological singularity, 
previous works [|l^ have shown that, in a suitable synchronous coordinate system (Gauss 
coordinates), the asymptotic behaviour of the string frame metric near the singularity reads 
(in any space-time dimension D) 

dsl ~ -dtl + E (-i5)'""(^n^f (^) , (5.15) 

a=l 

0(x, t) = 0(x, 0) + aix) \ogi-ts) , (5.16) 
where Ef[x) is some {D — l)-bein and where (1 < a < -D — 1) 




(5.17) 



In the Einstein frame this asymptotic behaviour reads 

dsl ~ -dtl + E(-^i^)'^"^"He- (x) dxy , (5.18) 



0(x, t) = <j){x, 0) + 7(x) log{-tE) , (5.19) 
where e'}{x) is proportional to Ef{x) and where 

D-l 

E Aa = 1 , 7 = ±v^ 



a=l 



\ 



l-EAJ. (5.20) 



a=l 



The "Kasner" exponents, in the string frame or Aq in the Einstein frame, can vary 
continuously along the singularity. The string parametrization is the most global one as it 
shows that aa runs freely over a unit sphere in R^^^ while the exponent for e'^ is a linear 
function of the "vector" aa- In the Einstein frame, the exponents Aq are restricted by the 
linear equality J2a^a = 1 and the quadratic constraint J2a^a < 1. A convenient geometric 
representation of these constraints, when D = 4, is to consider (in analogy with Mandelstam's 
variables s,t,u) that the 3 Aa's represent, in some Euclidean plane, the orthogonal distances 
of a point A from the three sides of an equilateral triangle (counted positively when A is 
inside the triangle). The quadratic constraint then means that A is restricted to stay inside 
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(1,0,0;0) 



(0,1,0;0) 




(0,0,1;0) 



+1 



■Y3 



(I o-^ 

(-l/3,2/3,2/3;0) 

FIG. 4. Geometric representation of the Einstein-frame Kasner exponents (Ai, A2, A3; 7) of Ein- 
stein-dilaton cosmological singularities. The three A 's, such that J2a Aa = 1 are the orthogonal 
distances to the three sides of an equilateral triangle. The constraint 7 = ib\/2 \/l — J2a 
stricts the representative A of Xa to stay on a two-sided disk circumscribed around the triangle. 
The pure-Einstein cases = 0) correspond to the circumscribing circle. In the spherically sym- 



metric case (X2 = X3), A runs over a bissectrix of the triangle. The basic parameter a of Eqs. (5.24) 
and ( 5.3i\ ) runs over a full real line (shown folded on the right of the figure ) and is mapped to the 
bissectrix via horizontal lines. 



the circle circumscribing the triangle. The sign ambiguity in Eq. ( p. 201 ) means that the 
parameter space is in fact a two-sided disk, namely the two faces of a coin circumscribed 
around the triangle. See Fig. ^. 

The link between the two seemingly very different parameter spaces (the string-frame 
sphere parametrized by the a's and the Einstein-frame two-sided disk parametrized by the 
A's) is geometrically very simple: one obtains the disk by a stereographic projection of the 
sphere, from the projection center (a^) = (1, 1, 1) (outside the sphere) onto the plane dual 
to this center (i.e. the plane spanning the circle of tangency to the sphere of the straight 
lines issued from a^). Algebraically (in any-space-time dimension D) 



aa{D -2) -a 



7 



2{D-2)(T 



(5.21) 



D-2-a ' ' D-2-a ' 

where a = (J^a'^a) — 1- Note that there is no sign ambiguity in the map — > A^. 
Coming back to the case D = 4, as aa runs over a unit sphere in R^, Xa covers twice the 
disk circumbscribed to an equilateral triangle. Each side of the disk is the image, via the 
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FIG. 5. Geometrical representation of the link between the two-sided disk of Fig. |^ on which 
the Einstein-frame Kasner exponents Xa live, and the sphere, (^a — ^> '^^ which the string-frame 
Kasner exponents {01,02,0^) live. The map Ua Xa is a stereographic projection, from the 
center (a^) = (1,1,1). The upper side of the disk comes from the projection of the polar cap 
a = iJ2a'^a) — 1 > located on the same side as a^, while the lower side of the disk comes from 
projection of the antipolar cap o" < 0. The points a+ on the polar cap, and a_ on the antipolar 
cap, are both projected on the same point A of the disk. The tangents to the sphere issued from 
cf touch the sphere along the circle limiting the disk (on which a = 0). 

stereographic projection, of a polar [a > 0) or antipolar [a < 0) "cap" of the sphere. This 
geometrical picture of the link between the a-sphere and the two-sided A-disk in represented 
in Fig. ^ It is useful to keep in mind this geometrical picture because it shows clearly that 
one can continuously pass from one side of the disk to the other, i.e. the sign of the Einstein- 
frame "Kasner" exponent of e'^ ~ (—tE)"'^^^ can change as one moves along the singularity. 
Physically, this means a very striking inhomogeneity near the singularity: the patches where 
7(0;) > (decreasing g'^ = e'^) will shrink (in string units) while only the patches where 
7(0;) < (growing g'^) can represent pre-Big Bangs. Finally, we note that Belinsky and 
Khalatnikov have proven that, in the generic case, the Einstein-frame Kasner exponents 
(Ai, A2, A3) must all be positive for the ultimate stable asymptotic approach to the singularity. 
This corresponds to A being inside the triangle of Fig. ^, i.e. in the hatched region of Fig. ^ 
However, this generic restriction does not apply in the special case of our toy model as we 
explain below. 
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C. Behaviour near the singularity in inhomogeneous, spherically symmetric solutions 



Let us now restrict ourselves to our toy model and study its cosmological-like behaviour 
near the singularity. First, we note that spherical symmetry will impose that two of the 
metric Kasner exponents (those corresponding to the 6 and ip directions) must be equal, say 
-^2 = -^3- Geometrically, this means that the Einstein- frame Kasner exponents run only over 
the intersection of the above two-sided disk with a bissectrix of the triangle. Second, we 
note that because of the vanishing, in spherical symmetry, of the crucial dreibein connection 
coefficient ei ■ (V x el), generically responsible for causing the expansion in the (radial) 
direction 1 to oscillate as t ^ 0^ ||3l|, a negative value of Ai, corresponding to an expansion 
in the radial direction, is allowed as ultimate asymptotic behaviour. As a consequence, even 
the portion of the bissectrix outside the triangle is allowed as asymptotic state. We shall 
prove directly this fact below by constructing a consistent expansion near the singularity for 
Ai of any sign. 

The mass-evolution equations (|3.5|) , (|3.6| ), together with Eqs. (|3.11| ), (|3.12| ) show that, 
starting from some regular event at the center, r = 0, where m = 0, the mass m{u,v) 
will grow if we follow a outgoing characteristic (i.e. a null geodesic) u = const. If this 
characteristic crosses the apparent horizon A, Eqs. ( |3.11|) , (|3.12| ) then show that m{u,v) 
continues to grow within the trapped region T{2m/r > 1) if one moves along either (future) 
outgoing or ingoing characteristics. This means that, within T, m{u,v) grows in all future 
directions. As the space-like singularity | B lies in the future of T, we see that the mass 
function will necessarily grow toward i3, and therefore, either will tend to a finite limit on 
i3, or it will tend monotonically to +oo as r — with fixed v (or u). Motivated by the 
Kasner-like behaviours Eqs. ( ^.18|) and ( 5.19| ) we expect the following generic asymptotic 
behaviour of m on the singularity (here written in (w, r) coordinates) 



^We consider here only the "non central" singularity, i.e. the part of the singularity which can be 
reached by future outgoing characteristics issued from the regular center r = 0. As shown by 
there exists also, at the "intersection" between the regular center and a singularity reachable 
only via ingoing characteristics. Note that r{u,v) tends to zero both at the center and at the 
singularity (central or non-central). 
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m(f , r) ~ Cmiv) — as r — * , f fixed (5.22) 



with some real f -dependent exponent a (v). [Here Iq denotes some convenient length scale.] 
It follows that yU = 2m/ r always tends to +oo on B (even when m tends to a finite limit). 
This suggests that we can consistently compute the structure of the fields near B by using an 
"anti-weak-field" approximation scheme in which /i ^ 1 (instead of the weak-field algorithm 



of section |IV A| where we used <^ 1). The expected domain of validity of such an "anti- 
weak-field" , or "cosmological-like" , expansion is sketched in Fig. ^. 
To leading order in this large- limit Eq. ( p.l5| ) yields 



1 f d(p \^ I d\ogm\ 
1\d\ogr)^ \d\ogr ) 

from which we obtain the following asymptotic behaviour oi (p on B 



(5.23) 



<j){r,v) ^2a{v)\ogr + C^{v). (5.24) 

We see now that the basic Kasner-like-exponent on B is half the coefficient of logr in 0(r, v) 
as r ^ 0, with fixed v. The Kasner-like exponent of m(f , r) is the square of that basic 
exponent. Moreover, the function Cm{v) appearing in the asymptotic behaviour of m(f , r) 
on the singularity is not independent from the functions a{v) and C^{y). Indeed, to leading 
order Eq. (|3.16| ) reads 

dlogm 1 dip d(j) 
dv 2 9 logr dv 

This relation is identically satisfied at order logr, and yields at order r° 



(5.25) 



dv dv 

Using now Eq. (|1|), written as 

f 9(3 \ 



aiogC^ dC^ 

a{v) — . (5.26) 



(5.27) 



y d log r J ^ \2 d log r J 
we further get the asymptotic behaviour of l3{v, r) 

p{v,r) a^{v)logr + Cfi{v) . (5.28) 

Finally, it is easy to see that the asymptotic behaviour of the metric on B in {v, r) coordinates 
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^^2^ 2a>)^2C,M ('M"^''^rf^2^2e^,w (Ly^''\ydr + r'du;\ (5.29) 



is indeed of the expected Kasner form (|5.18|) by introducing the cosmological time 



(-0 - '•(-) . (5.30) 



The transformed metric reads, to dominant order as r — 0, i.e. t , 



2C„ 



0~7(t;)log(-t), (5.32) 

The basic parameter a{v) runs over the real hne — oo < a < +oo. One easily checks that 
the relations ( |5.2(j| ) are satisfied. Some cases having a particular significance are illustrated 



in Fig. ^: When a = 0, one gets (Ai,A2,A3) = (— and 7 = which describes 
a Schwarzschild-like singularity (when viewed with a cosmological eye). When a = ±1, 
one gets (Ai,A2,A3) = (o, |, |) and 7 = ±1 which is locally similar to the scale-invariant 
solution |2^ also described earlier. When a = ±-\/3, one gets (Ai,A2,A3) = (|) I? l) and 
7 = ±2/-\/3 which describes an isotropic cosmological singularity of the type of the exact 
solution [|30|, |]18[ we discussed above. As illustrated in Fig. H, when a runs over M its 
image, defined by Eq. ( |5.33|) , runs twice over the bissectrix of the Mandelstam triangle 



(intersected with the circumbscribed disk). The only pure-Einstein case (7 = 0) reached at 
finite a is the Schwarzschild-like case a = 0. In principle, another pure-Einstein behaviour 
occurs when |a| — > 00, leading to (Ai, A2, A3) = (1,0,0) i.e. ds"^ ~ —dt"^ + dp^ + doj"^ = 
—dr"^ + d^"^ + d9^ + sin^ 9 dip^. This would be an Einstein-cylinder-like universe. We believe, 
however, that \oi{v)\ will stay bounded (|a('y)| < ttM < +00) as v runs over B and that the 
only Einstein-like case will be crossed when a{v) changes along B. 

We have confirmed our asymptotic analysis of the structure of the singularity (based on 
the Ansatze (|5.22| ), ( ^.24| ) and ( |5.28| ) in two ways. First, we have verified that it is consistent 



with the partial, but rigorous, results given in Ref. [^. Namely, in {u,v) coordinates, our 



34 



Ansatze are compatible with the result that the function r^{u, v) is on B if we posit that 
the conformal factor behaves like Vl? ~ fO^i-")-^ (while m ~ r~°^*^^^). A stronger check of our 
assumptions is obtained by going beyond the leading order approximation, thus extending 
some results of [0. In analogy to section P^V A| , where we set up a complete perturbation 
algorithm in the weak- field domain, we have shown that, starting from the leading-order 
terms Eqs. (|5.22| ), (|5.24|) and (|5.28| ), containing 3 arbitrary "seed functions", C0(f), 
Cpiy) [from which Cmiy) can be determined using Eq. ( p.l6[ )], it is possible to set-up an 
all-order iterative scheme generating a formal solution of the spherically symmetric Einstein- 
dilaton system in the strong- field, cosmo logical-like domain. Note that, Cj3{v) being a pure 
gauge function (it suffices to introduce v' = J e'"'^^'"^ dv to gauge C^(f) away), the formal 
solution thus generated contains two physically arbitrary functions of v, which is indeed the 
freedom that should be present in a generic solution Details of this scheme are given in 
Appendix 0. We only mention here the general structure of the expansions so generated 
near the singularity, i.e. as r — > 0, with fixed v. (see Eq. (|A34|) in the appendix) 



^(^y^r) = Mv,r)+ a?ir"("'(^)+i)r2-p/<^)(logr), (5.34) 

n,m,l 

(3{v,r) = Po{v,r) + ^ aj^i r"("'^")+'^ r^^ Pl''\logr) , (5.35) 



m{v,r) = mo{v,r) + ^ r"(" r''" Pr^(logr) . (5.36) 

n,in,l 

Here n, m, / are integers with I < m, a'nmi ^"^^ coefficients of the polynomials P/'^'' are 
functions of f, and (0o, /^o, ""^o) are the leading-order terms given by Eqs. (|5.22| ), ( ^.24| ) and 
( ^.28| ). This expansion actually contains two intertwined series: an expansion in powers of 



X 



and a more complicated series in r^"^ (logr)'. The second expansion is linked to 
the w-gradients of the seed functions a{y), Cfj,{y), Cp{v), while the first expansion is present 
even in the simple "homogeneous" case where a, and C/3 do not depend on v. The fact, 
exhibited by Eqs. (|]34D-(|]36D, that the expansion in the homogeneous case proceeds along 
powers of r"^+^ confirms a recent conjecture by Burko [p5|] . 



^'^Here, contrary to what happened above, the variable v varies only on a half-line, +00] where 
v^, marks the birth of B. As was said above two functions on a half-line (or one function on the full 
line) is the correct genericity in our toy model. 
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D. Exact homogeneous cosmological solutions 



Let us consider in more detail the "homogeneous" case where the seed functions a, and 
Cj3 have no spatial variation along the singularity. As it turns out, it is possible to resum ex- 
actly all the terms of the "homogeneous" expansion. Indeed, the analog of the Schwarzschild 
solution (t- independent, spherically symmetric solution) for the Einstein-dilaton system has 
been worked out analytically long ago by Just (see also using the special gauge 

where googrr = — 1- Contrary to the case of the Schwarzschild solution, this solution cannot 
be continuously extended (through a regular horizon) down to a cosmological singularity at 
r = 0. However, it is easy to see that the following cosmo logical-like background (which is 
related to Just's original solution by formally extending the radial variable r = Ix "below" 
the curvature singularity at r = a) is still a solution of the Einstein-dilaton system: 



- ( dx' + x'+' (1 - x)'-' 

X J W — X J 

1 — X 



(5.37) 



X 



±Vl^ log . (5.38) 



Here the (formerly radial) variable x varies between and 1 and is "time-like", while the 
(formerly time) variable t is "space-like" . The cosmological universe evolves from a Big Bang 
at X = to a Big Crunch at x = 1. It has two arbitrary parameters, a scale parameter /, 
and the dimensionless 6, — 1 < 6 < 1. Near the singularity at x = 0, the link between h and 
the parametrization used above is 



Note that the behaviour near the other singularity at x = 1 is obtained by changing the 
sign of h and by changing a into —1/a. Some cases of this homogeneous solution are of 
special significance: 6 = 1 is Schwarzschild (with x = r/2m), 6 = belongs to the class of 



scale- invariant solutions discussed above and h = —1/2 interpolates between a locally 
isotropic cosmological solution (a = =FV^) at x = and an anisotropic one (a = ±l/v^) 
at X = 1. Note that, in spite of the homogeneity and isotropy at x = 0, the latter special 
solution differs from the (everywhere) homogeneous-isotropic solution described earlier. 
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E. Discussion 



Ideally, the two formal expansions we have constructed, the weak-field one Eqs. ( |4.1| )- 
([4 .31), and the strong-field one Eqs. ( |5.34| )-( p.36|) , should match at some intermediate hyper- 
surface, like the apparent horizon 2m /r = 1 which looks like a natural borderline between 
the two near domains. If this matching were analytically doable, it would determine all the 
(so-far) arbitrary "seed functions" of the strong-field scheme in terms of the unique arbitrary 
function of the weak- field one, namely the asymptotic waveform f{v) (or, equivalently, the 
news N{v) = f'{v)). However, it is clearly too naive to expect to perform this matching 
perturbatively: neither of the two expansions is expected to be convergent |^ (they are prob- 
ably only asymptotic). Even if they converge on some domain they probably both break 
down before reaching a possible overlap region where they might be matched. At this stage, 
we can only state that, in principle, all the seed functions of the strong-field scheme are 
some complicated non linear and non local functionals of N{v). It would be particularly 
interesting to study the functional dependence of the Kasner exponent a{v) on N{v). By 
causality (i.e. a domain of dependence argument) we know that <y{v), at advanced time 
V, depends only on N{v') on the interval v' < v. When starting from a generic N{v) "of 
order unity", we expect that the resulting a{v) will also be of order unity. A physically very 
important issue is the sign of a{v), i.e. the sign of 7(f). Indeed a > means a decreasing 
0, while a < means that (p grows near B. Let us note that, to lowest-order of weak- field 
perturbation theory, the value of 0(f , r) at the (regular) center (r = 0) is 

f{v)-f{v-2ry 



lim 4>i(v, r) = lim 



+ 



<Pi + 2f'{v) = <Pi + 2N{v) . (5.40) 



Here (pi denotes the background value at past infinity. From this result we expect that, 
if N{v) is a simple Gaussian-like wave packet, N{v) = Aexp [— (t> — t>j)^/£^], with a large 
enough dimensionless amplitude for leading to collapse, i.e. \A\ ^ 1, the local values of 
near the collapsing central region will, at first, grow if A > 0, and decrease if A < 0. In 
this simple case, we therefore expect a{v) to have the opposite sign of A near the "central" 



^^The possible presence of a finite number of BKL oscillations near B suggests that the formal 
strong- field expansion has very bad convergence properties. 
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region of B, i.e. for v ^ Vi where the gravitational instabihty sets in. But, further away 
on B, i.e. for v ^ Vi, the sign of a{v) can change. We also expect that, as f +00, 
a{v) will tend to zero, corresponding to a Schwarzschild-like singularity (in the case of well 
localized incoming N{v) packets). This conjectured link between the sign of a{v) and the 
sign of N{v) is confirmed by the exact isotropic solution Eq. (|5.1| ) (with a positive growing 
N{v) = v^/f^, and a = — v^), as well as by the scale invariant solution Eq. ( [4.40|) . Some 
numerical calculations seem to confirm the general picture we propose (in particular the 
interesting possibility that the sign of a{v) changes several times on as f varies, before 
reaching a Schwarzschild-like asymptotic regime |a(f)| — as f — +00). We plan to study 
in more detail these issues in a future publication. 



VI. A BAYESIAN LOOK AT PRE-BIG BANG'S "FINE-TUNING 



From its inception |^ it was pointed out that a successful pre-Big Bang scenario must 
rely on a "reservoir" of inflationary e-folds during its perturbative phase. This is given by two 
small numbers, the initial curvature scale in string units, Hiig, and the initial string coupling 
constant gf = e"^^ 1. [In this section, the index i is used for labelling "initial" quantities.] 
The need of large (or small) numbers has been recently discussed at length [|TB|, ||TP|] and 
used to criticize the naturalness of the PBB scenario. In particular, it was pointed out 
T8| , |jT9| that, as soon as one goes beyond the simple spatially-fiat, homogeneous cosmology 



framework, the total duration of the perturbative dilaton-driven phase is finite so that the 
resolution of the homogeneity /flatness problems requires 



9i 



<A./2< 10-26, t^>t^>iQ^\ (6.1) 



Here, Lj denotes the spatial homogeneity scale of the PBB universe at the beginning of its 
inflationary phase, and Hi ^ L^^ its time-curvature scale (Hubble parameter). 

In this Section we shall systematically work with string-frame quantities, even when 
we refer to results discussed in previous sections in the Einstein frame. For instance, the 
time-scale ii, characterizing the rate of variation of the news functions around the advanced 
time Vi and leading to gravitational instability, is now (locally) measured in string units. In 
any case, we are essentially working with dimensionless ratios which are unit-independent: 
{^h/^'s)s = {^h/^s)e = dii^i/ ^pia.nck)E- We wish to emphasize here that, by combining our 
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stochastic-like instability picture with a Bayesian approach to the a posteriori probability of 
being in a position of asking fine-tuning questions, the issue of the naturalness of the PBB 
scenario is drastically changed. By "Bayesian approach" we mean taking into account the 
selection effect that fine-tuning questions presuppose the existence of a scientific civilization. 



As emphasized long ago by Dicke |39| and Carter [40|, civilization-related selection effects 
can completely change the significance of large numbers or of apparent coincidences. Linde 
and collaborators explored several aspects of the Dicke-Carter "anthropic principle" within 
the inflationary paradigm, and emphasized the necessity to weigh a posteriori probabilities 
by the physical volume of inflationary patches [^. We shall here follow Vilenkin ||4^, and 



his "principle of mediocrity" , according to which the unnormalized a posteriori probability of 
a random scientific civilization to observe any values of the PBB parameters Qi, Li and Hi is 
obtained by multiplying the corresponding a priori probability by the number of civilizations 
(over the whole of space and time) associated with the values gi, Li and Hi. AfdviOi, Li, Hi). 

A. A priori and a posteriori probability distributions for gi and Hi 

This approach can be applied to our case if we think of the initial past-trivial string 
vacuum as made of a more or less stochastic superposition of incoming waves (described by 
complicated news functions N{v, 9, (p), A^+(f , 9, ip), A^x {v, 9, ip) having many bumps, troughs 
and ramps). This stochastic bath of incoming waves will generate a rough sea of dilatonic 
and gravitational fields. If the input dimensionless wave forms can reach values of order 
unity, we expect that the local conditions for gravitational instability will be satisfied at 
several places in space and time. This will give birth to an ensemble of bubbling baby 
universes, with a more or less random distribution of initial parameters Qi and Hi, and with 
initial spatial homogeneity scales Li > H^^. Indeed, the analysis of the previous sections 
has shown (in our toy model) that gravitational instability will set in on a spatio-temporal 
scale ~ ii when a rising wave of news function grows by 0{1) on an advanced-time scale 



ii. From our variance criterion Eq. ( 4.25 ), and the rough validity of weak- field perturbation 
theory nearly until the sharp transition to a cosmological-type behaviour, the work of the 
previous sections has shown that the initial, advanced-time scale ii is propagated via ingoing 
characteristics with little deformation down to the strong- field domain (see Fig. ^), where it 
appears as a spatial homogeneity scale, i.e. 
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(6.2) 



From the leading order result ( |5.40| ) we expect the local value of the Hubble parameter in 



the corresponding cosmological-like bubble to be 

V^~|iV'(t;,)|~M!^ (6.3) 

where Vi denotes the advanced time at which the instability sets in. Combining ( |6.2|) and 



( p73D we get 

HiU^\N{vi)\>l. (6.4) 

These rough formulae show how, in principle, given the stochastic properties of the dimen- 
sionless news functions, one could deduce the distribution of Hi and Lj, naturally constrained 
hy Hi Li ^ 1. The corresponding distribution of Qi is a priori independent from that of 
ii ~ Lj, being linked to the presence of a slowly-varying ("DC") component Nq in N{v) (by 
contrast to the local variations of N{v) leading to instability): N{v) = A^o + fluctuations. 
The DC component A^o corresponds to ramps in /(f), f{y) = Nqv + fo + fluctuations, i.e. 
to shifts of (pi : (pi (pi + 2Nq. Finally, we can consider that the initial distribution of the 
news functions deflnes an a priori probability distribution for Qi, and Hi, 

dpprior = Wi{gi, Li) dgi dLi , (6.5) 

with Hi generically of order L~^ (because |A^(fj)| ~ 1 is the threshold for instability). 

An important remark must be made here. The two basic parameters we are talking about. 



say Qi and Li Hi , precisely correspond to the two global symmetries, Eqs. (|2.9|)-(|2Tl|) 



and Eq. (|2.12|) , of the classical string vacua: a constant shift in 0, and a global coordinate 



rescaling. If we consider that the initial state of string theory is classical (rather than 
quantum), these symmetries mean that no particular values for (pi or Li are a priori preferred. 
One might even expect a "flat" distribution compatible with these global symmetries. 



^ Hassicalflat ^ ^ ^ ^ ^ ^ . _ (g,g) 

gi L 



Evidently, the problem with such a flat distribution is that it is non-normalizable. Some cut- 
offs are needed to make sense of such a flat prior distribution but, while there are natural 
strong coupling/curvature cut-offs (the limits of validity of our approximation), it is not 
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so easy to find natural small-curvature/coupling cut-offs. If we appeal to string theory for 
providing cut-offs for gi and Lj, arguably the prior distribution ( |b.5D should be determined by 
the conjectured basic symmetries of string theory, i.e. by S and T-duality. ^-duality suggests 
that one should work with the complex quantity z = ^ + ^ p- and require modular invariance 
m z = X + iy. This selects the Poincare metric ds^ = {dx^ + dy^^jy^ as being special, and 
thereby defines a preferred measure in the fundamental domain (key-hole region) of the 
(x, plane: rf/i = ^dxdy = dxdy/y"^ = dx\dy^^\. Integrating out the angular variable 
6 over < 6 < 2tt (which is the correct key-hole range when g"^ < Att, i.e. y > 1) we end 
up with a preferred probability distribution for g : dp{g) oc dg"^ oc g dg, considered in the 
weak-coupling region g 1. A similar argument based on T-duality selects as preferred 
probability distribution for the spatial length scale L : dp{L) oc dL'"^ oc L^^ dL where 
L = L/is is considered in the weak-cr-model-coupling region L > 1. In conclusion, this type 
of argument would suggest the following normalizable, factorized prior distribution 

c^P^OT^'^"''' OC (gi dgi) L^^ dU oc gi dgi Hi dHi , (6.7) 

with < gi< 1,0 < Hi = Hi/M, = Lr^ < 1. 

For the sake of generality, we wish to leave open the nature of the prior distribution, and 
thus consider the general class of prior probability distributions: 

gt dgi H'^ dHi , (6.8) 

with arbitrary powers a and b. 



Any given initial distribution, say (|6.8|) , will generate a corresponding ensemble of PBB 
infiationary bubbles. We assume, as usual, the existence of a successful "exit" mechanism 
by which dilaton-driven infiation, with growing g and H finally "exits" into a standard hot 
Big Bang, i.e. a radiation-dominated Friedmann universe. Current ideas about how this 
might happen assume that the Hubble parameter reaches values of order of the string 
mass Ms, i.e. Hg = H^/Ms ~ 1, before the string coupling constant reaches values of order 
unity. If the opposite were true ((? ~ 1 occurring before ~ 1), it is likely that quantum 
fiuctuations would take over before H reaches the string scale, causing inhomogeneities to 
grow so large that the infiationary process "aborts" before a baby universe is born. In any 
event, we shall discard this possibility, assuming that the resulting cosmo logical universe 
would not evolve into anything able to harbour scientific civilizations. 
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After H reaches 1 at g <^1, the evolving universe is expected to be entering a so-called 
stringy De Sitter phase, i.e. a phase during which and H are constant and given by some 



fixed-point values of the order of M^: |0| = xM^, and \H\ = y Mg with a; ~ ?/ ~ 1. [See 
for some suggestions on how to implement this mechanism.] Only the ratio between x and y 
enters our present phenomenological discussion. In keeping with a notation used in previous 
papers on PBB phenomenology, we introduce the (positive) parameter /3 by 

|(ln^)-| = ^101 = /? if = const. (6.9) 

One finally assumes that this stringy De Sitter-like phase ends when g = e'^^'^ reaches 
values of order unity, causing the amplified vacuum fiuctuations to reach the critical density. 
After this moment, g should remain fixed at about its present value, either because of a 
nonperturbative 0-dependent potential, or possibly because of the attractor mechanism of 
Ref. 0. The fixing of g marks the end of the stringy modifications to Einstein's theory, 
and the beginning of standard Friedmann-like cosmology. As we assume that g is attracted 
towards a unique fixed point, the resulting Friedmann universe will have no free parameter 
which could make it different from our universe, except for its initial homogeneity scale Lf, 
where the index / refers to the final state of the stringy phase. Therefore, the number of 
civilizations in the resulting Friedmann universe will be proportional to its total volume, at 
least for the large enough (in fact, old enough) universes that can harbour life 

KU9i, H,) CX V/((?„ H,) e[Vf - Knin] . (6.10) 



One assumes here, a la Dicke |3^, that the time span during which civilizations can occur, 
being constrained by the lifetime of stars, is fixed. The minimum scale Vmin for a life- 
harbouring universe is not known precisely, though it corresponds probably to universes 
whose total lifetime before recollapse is a few billion years. In conclusion, the a posteriori 
probability for a random scientific civilization to ask fine-tuning questions about the values 
of gi and Hi is 

dppost oc Afcivigi, Hi) dpprioT oc Vf{gi, Hi) 9{Vf - Vmin) dpprioT ■ (6-11) 

B. Computation of the a posteriori probability distribution in 4 dimensions 

The computation of the final volume V/ can be separated as 
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Vf = L} = Lfy^^zl (6.12) 

Here, the index i refers to the (cosmological) time ti of beginning of dilaton-driven inflation, 
the index 1 refers to the time ti of transition between dilaton-driven superinflation and 
a De Sitter-hke phase, and Zs = exp {H (t2 — ti)) denotes the total expansion during the 
De Sitter phase. The ending t2 of the De Sitter-like phase is assumed to take place whence 
g2 = git2)r^Oil). 

Using the leading power-law behaviour of the cosmological evolution during the dilaton- 
driven phase, Eq. ( |5.15|) , and imposing that this phase ends when Hi ~ (— ti)^^ ~ Mg, one 



gets 



where Sa = ai + ^2 + 0:3 is the sum of the string-frame Kasner exponents. The a's are 
negative (expansion) and are constrained by Eq. (|5.17| ) . In the result (|6.13|) we have neglected 



a factor of order unity linked to the fact that the dilaton-driven inflation is generically 
anisotropic (ai 7^ 02 7^ cts) and therefore that the various Hubble expansion rates will 
reach the string scale at slightly different times. We assume here that a nontrivial basin of 
attraction from anisotropic expansion toward an isotropic stringy De Sitter-like phase exists 



as it is the case in examples |43 



During the subsequent stringy De Sitter phase, the growth of g is related, according to 
Eq. (|6.9[), to the spatial expansion hj g oc so that 



«i \9i, 



' ^g'lK (6.14) 



where we used (72 ~ 1- On the other hand, gi = g(ti) is expressible, through the use of 
Eq. ( I5.16D , giving the evolution of g during the dilaton phase, as 



(6.15) 



\gij \U 

Finally, we express, using Lj ~ Hj^^, the volume at the beginning of standard cosmology in 
the form 

Vf{g,, H,) ~ H-^ [gi{g,, %)]-^/^ e{l - g^) , (6.16) 
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with 



g,{g,,H,)^g,Hr^, (6.17) 
and where the exponents A and B are positive, and are given in terms of the a's by 

A = 3-Sa, 5 = ^-isa. (6.18) 



If we assume, for illustration, an a priori distribution of values of gi and Hi of the form (|6.8|) , 
we get for the a posteriori probability distribution of gi and Hi 

dp^ost oc gt Ht^ [g^ H^^^^^ 6(1 - g,{g,, Hi)) e{Vf - V,nin) dg, d% , (6.19) 

which can be written as 

dp^ost oc glH;^ Oil - g^) Oiyj - V.,i„) dg^ d% . (6.20) 

In the last form we have expressed the measure dp^o^t in terms of the two independent 
variables gi and Hi (this is convenient because of the constraint gi < 1 carried by the step 
function). The exponents a and A appearing there are 
3 

a = a — — , (6-21) 

/-^ 

— 11 

A = A- {a+l)B -h = ?> + -{a-l) (Sa) - - (a + 1) - 6 . (6.22) 

The numerical values of the exponents a and A play a crucial role in determining the a 
posteriori plausibility of our universe having evolved from the seemingly "unnatural" values 
( |6.1|) . Indeed, if A > 1 the posterior probability distribution for Hi is peaked, in a non 
integrable manner, at Hi = 0. Therefore, if A > 1, it is natural to expect (within the 
present scenario) that the initial homogeneous patch, whose gravitational instability led 
to our universe, be extremely large compared to the string scale [and that the value of 
gi = H^ gi with B > and gi < 1 he correspondingly small] . Most scientific civilizations 
are bound to evolve in such universes. 

Note that the duality-suggested values a = 1 and 6 = 1 of Eq. (|6.8|) lead to A = 1 
(independently of the a's) for which such an a posteriori explanaiion works. From this point 
of view one can argue that, within string theory, it is "natural" to observe very small numbers 
such as in Eq. (|6.1|) . The exponent a also plays an important role. Indeed, if a > —1, i.e. if 
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/5 > ^ , (6.23) 

a + 1 

the a posteriori probability distribution for gi is integrable over its entire possible range 
(0, 1). This means that most scientific civilizations are expected to observe values gi = 0{1), 
corresponding to a small number of e-folds during the stringy phase: ~ gi^^'^ = This 
is a phenomenologically interesting case, because it means that various interesting physical 
phenomena taking place during the dilaton-driven phase (such as the quantum amplification 
of various fields) might leave observable imprints on cosmologically relevant scales. [A very 
long string phase would essentially iron out all signals coming from the dilaton phase.] It is 
interesting to note that, in the string-duality-inspired case ( |6.7|) , the great divide between a 
short stringy phase (a > —1), and a very long one (a < —1, implying that gi is peaked in 
a non integrable manner at gi = 0), lies at f3 = 3/2, which played already a special role in 
previous phenomenological studies 0. 

Though the case A > 1, a > —1 appears as the conceptually most interesting case for 
the pre-Big Bang scenario, other civilization-related selection effects might also render the 
case A < 1 viable. In this case, the factors H~^dHi in Eq. (|6.2CI|) suggests that Hi should 
be of order unity (i.e. Li ~ H^^ ~ £s). However, in such a case, one must take into account 
the factor 6'(V/ — Vmin) (which could be neglected in the previous discussion). This factor 
means that scientific civilizations should a posteriori expect to find themselves in the smallest 
possible universe compatible with their appearance. This probably means that they should 
also expect to see inhomogeneities comparable to the Hubble scale, and to have appeared 
very late in the cosmological evolution, just before the universe recollapses. This seems to 
conflict with our observations. In conclusion, within this model, and limiting our discussion 
for simplicity to the class (|6.8| ) with a and h integers, one finds the favorable situation A>1 
realized (i.e. no a posteriori unnaturalness in observing very small gi and Hi) when either 
a = and 6 < 2, or a = 1 and 6 < 1. However, if a > 2, i.e. if the a priori distribution 
function of gi vanishes faster than gf for gfj — 0, the pre-Big Bang scenario has to face a 
naturalness problem. [This is the case for quantum fluctuations of gi, as discussed at the 
end of the next subsection.! 
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C. Computation of the a posteriori probability distribution in D dimensions 



Assuming that the results discussed so far are quahtatively correct in higher dimensions, 
we have generahzed the above considerations to the case in which the initial string vacuum 
is extended to d spatial dimensions (say = 9), and where, through gravitational instability, 
3 of these dimensions are expanding and 6 are collapsing down to the string scale. This 
leads to introducing 4 special times: ti (beginning of dilaton phase), ti ~ —ts (end of 
the dilaton-driven power-law evolution of space), t2 (when the radius of curvature of the 
collapsing dimensions become 0{is))i and (when = 0{1)). In this scenario there are 
two stringy De Sitter-like phases: a first one during which the 3 "external" dimensions 
grow exponentially a oc exp(ifit), while the {d — 3) "internal" ones shrink exponentially 
a oc exp {—Hit) (and while = 2[3iHi = 2(3-1 Hi), a second phase during which the 
{d — 3) string-scale-curved internal dimensions have frozen while the 3 external ones continue 
to grow exponentially a oc exp {H2 1) (and while = 2/^2 H2). This scenario gives a result of 
the same form as above, namely 



Vf{g„Hi)r^ Hr^ g^ig^H, 



'3//32 



[^-92] 



(6.24) 



with 



92(91, Hi) = giH- 



-B 



(6.25) 



and where 



A = d - 3a - (d - 3)a + (1 - a) (3i I ^ - 'L^ 



(6.26) 



B 



■ a 



{d-3)a + (l-a) pi. 



(6.27) 



Here one has assumed for simplicity that, during the dilaton phase, the 3 external dimensions 
grow with the same exponent a < 0, i.e. proportionally to (— 1)°, while the {d — 3) internal 



12 



This string-frame Kasner exponent a should not be confused with the basic parameter a{v) 



introduced in Sec. V C 
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ones decrease with the same exponent a > 0, proportionally to (— Then, one can derive 
the analog of Eq. ( |6.2(J| ) with gi replaced by g2 and with exponents 

a = a-|- ~A = A-{a + l)B~h. (6.28) 

Once more, A>1 leads to no fine-tuning and a > —1 to a short string phase. We shall not 
attempt a complete phenomenological discussion of the range of values of a and h compatible 
with A > 1. (This would imply exploring the parameter space of allowed values for the 
Kasner exponents a and a.) Let us only note here that the introduction of more spatial 
dimensions generically helps because A (and A) now contains the contribution +d = +9, 
say, instead of the previous +3. In this case, even distribution functions vanishing faster 
than gf (for g^ 0) can make very small values of gi be a posteriori preferred. The 
power-law enhancement brought by the volume factor V/ would be inadequate, however, 
for compensating an exponential suppression as g^ 0, say oc exp {—c/gf). In particular, 
the a priori distribution of gi and Hi cannot be the one expected (from the Lagrangian 
C ~ e~'^(V(py ~ (V(7~^)^) for quantum fiuctuations of g on scales Lj ~ around 
the trivial ground state. Indeed, this distribution vanishes for small gi and small Hi as 
exp {—cgl"^ H^"^) with c of order unity. This confirms the standard idea of the pre-Big Bang 
scenario, according to which the initial state should be a classical string vacuum, i.e. an 
arbitrary classical solution of the low-energy field equations. 



VII. CONCLUSION 



We have proposed a new, concrete picture for the typical initial state of the Universe 
which can give rise to a cosmology of the so-called pre-big bang type. We conceived this 
state as a very generic, past-trivial string vacuum, i.e., classically, as the most general so- 
lution of the low-energy tree-level string action endowed with a perturbative region in the 
asymptotic past. Such states can be parametrized (in three spatial dimensions) by three 
dimensionless "news functions" : a scalar news N{v, 6, (f) and two helicity-2 news N^{v, 6, ip) 
and A^x {v., 0, (p). By studying the simple spherically symmetric tensor-scalar model we have 
argued that such a generic "in state" will be gravitationally unstable if (any of) the news 
functions varie(s) by something of order unity on some scale Af ~ ^. More precisely, in 
the spherically symmetric model, we found, at lowest order of perturbation theory, that the 
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criterion for gravitational instability is given by the supremum over Vi and f 2 of the variance 
of the function N{x) over the interval [^1,^2], see Eq. ( [4.25[ ). When this supremum exceeds 
a threshold of order unity, the wave packet "sent" between the retarded times vi and f 2 will 
become gravitationally unstable when, later, it becomes focussed at the center We have 
verified that the infinite region between past infinity and such a center can be described 
by weak-field perturbation theory. Then, near the center, in a region of spatial extension 
Li ~ Av = V2 — Vi, one can abruptly shift from a weak-field description to a cosmologi- 
cal one. One can then describe the further evolution of the patch Li by a complementary 
cosmological-like expansion of the type introduced by Belinsky, Khalatnikov and Lifshitz 
(BKL) and recently developed, for the PBB scenario, in |]12|. We have proven to all or- 
ders of iteration (in our toy model) that the presence of the dilaton ensured the consistency 
of a (non oscillatory) BKL expansion in log-corrected powers of the distance to the future 
(space-like) singularity. This result confirms the smoothing behaviour of a dilaton-driven 
superinflationary evolution. In the string frame, gravitationally unstable patches of vac- 
uum, with homogeneity scale Lj, will expand (if > 0) into homogeneous patches on much 
larger scales. These blistering patches are surrounded by contracting regions where < 0. 
Many such dilaton-driven pre-Big Bang inflationary bubbles, looking like closed universes, 
can blister off a generic past-trivial string vacuum. By taking into account the selection 
effects linked to the presence of scientific civilizations, we have shown that it can appear a 
posteriori probable to observe that our universe comes from a very large fluctuation Li ^ is, 
with a very small Qi. In particular, we pointed out that string-dualities suggest an a priori 
distribution function for Qi and Lj which is phenomenologically encouraging in predicting a 
short, intermediate De Sitter-like stringy phase, and the a posteriori natmalness of observing 
Li ^ is and gi <^ 1. 

Our scenario contains new features that might suggest new observable signatures of the 
pre-Big Bang idea. In particular, we expect the dilaton-driven phase to be generically 



In absence of spherical symmetry, there may exist many "local centers" where gravitational 
instability sets in. They should each correspond to an annular section of X~ with a 0(1) fluctuation 
of the news looking roughly spherically symmetric when viewed from some center and in some 
suitably boosted Lorentz frame. 
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anisotropic. It will be interesting to study the possible fossils of this anisotropy, especially 
in the presence of a short string phase. Our work suggests also new types of criteria for 
the gravitational instability of the Einstein-plus-scalar system. Concentrating on some non 
cumulative measure of the variation of the news might lead to new theorems of the type of 



those of Ref. [H, HI. 

One of the most interesting well-posed technical problems suggested by our work is the 
study of the link between the incoming news function and the spatial dependence (along the 
singularity) of the Kasner exponents. We hope to study this problem in the near future. 
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APPENDIX A: 



In this appendix we will establish the structure of the perturbative expansion near the 
singularity in the spherically symmetric model analysed in the text. We found it more 



convenient to use the (f , r) coordinate system introduced in Sec. |T|. The field equations 
read 



dvdr r \dv , 



+ e 



dr 



0, 



(Al) 



^ or J 



r I 0(p 
4 



dr 



2m\ / d(f)^ 
r J 4 \ dr , 



(A2) 



Let us expand the different fields {(f), (3, m) around the leading-order solutions (|5.22|) , (|5.24| ) 
and (lOSi ) 



0o(w,r) = C^{v) + 2a{v) logr , 
Po{v,r) = C/siv) + a^{v) logr , 
mo(t;,r)=C^(t;)r-"'(-), 



(A3) 
(A4) 
(A5) 
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and define the corrections to it 0, (3 and m by 



0(t;,r) = (f)o{v,r) + (piv,r) , 
m{v, r) = mo{v, r) (1 + rn{v, r)) . 



(A6) 
(A7) 
(A8) 



Introducing the variable = r, and denoting differentation with respect to f by a prime 
Eqs. (|Al|)-(^ take the form: 



S^{v, A) + L^{v, A) + NL^{v, A) , 

' V(l + A)e2^ + 2C;e2^ + 2ae^''e("'+i)" 



2a 

e 



m 



1 

2Cn 



,2A 



(a2+l)A 



'^2 



aA 



+ 



9A2 aA 



2a e("'+^) ^ (e^ - 1 - /5) - 2C„ m 



92 



aA -"Ua 



5« = = 



= Sfsiv, A) + A) + NL^{v, A) , 



(A9) 
(AlO) 

(All) 



(A12) 

(A13) 
(A14) 



Sry 



NLr 



2„(q2+i)a 



a e 



2CfYi 

„(a2+i)A 



, 2 



a e 



(a2+l)A 



2a 



m 
2 



8C„ 



dx -4 Ua -""^ Ua 



T IsA 



(A15) 
(A16) 
(A17) 



In the R.H.S.'s of Eqs. (|A9|) , ( |A13|) and ( |A15|) we have indicated by S the field-independent 
terms ("source terms"), by L the terms that are linear in the fields, but subleading with 
respect to the L.H.S. in the limit r — > (A ^ —00), while by NL we refer to generic 
non-linear contributions. 

By examining the source terms and the subleading linear terms, we see that they all 
contain either an explicit factor e^"^^^-''^ = ^(q^+i) q^. ^ factor e^^ = r^. To keep track of the 
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presence of these factors (which tend to zero at the singularity) let us introduce two formal 
bookkeeping parameters: (Xi associated with each occurrence of e^"^"*"^)^ in 5* and L and, 
likewise, o"2 associated with e^'^. In the end we shall set ai = (J2 = 1, but the introduction of 
these parameters will be useful to delineate the structure of the expansion of the fields near 
r = 0. 

Introducing the notation if) = (0, /3, m), considered as a column matrix, and denoting Dx 
the following derivative operator with respect to A, 



^ di ^ 



(A18) 



-a{v)dx dx 
\^ +a{v)dx dx J 

we can rewrite Eqs. ([A9|), ( |A13|) and ( |A15| ) in the compact matrix form: 

Dx^P = ai + (T2 + ^1 + ^2 if + NL^ . (A19) 
The nonlinear terms in ( |A19|) depend linearly on cxi and NL^{il),dx'4'id'x4^i<^ii<^2) = 



NL^ + ai NL^l^ + (J2 Nlf. The differential matrix system ( [AT9| ) can be rewritten as an 



integro-differential system by introducing the "retarded" Green function G(A, A'), i.e. the 
unique inverse of Dx which vanishes when A' > A 

^ = G * (ai + a2 + ^1 L^i^ + ^2 + ^^^) ' (^20) 



where we indicate by * the operation of integrating over A'. This choice of boundary con- 
ditions is imposed by consistency: We defined the background solution ( |A3|) -(|X5|) as giving 
the asymptotic behaviour of (0, (3, m) when r — 0, i.e. A —00 . We must therefore require 
that the deviation field ipiX) vanish when A —>■ —00. 

The system ( |A2CI| ) can be formally solved by successive iterations, leading to a solution 
ijj in the form of a double power series in the bookkeeping parameters o"i and cr2: 

^= E <<^nn^. (A21) 

n,m>l 

For instance, the first order (in ai and (72) corrections V'l = (^I'^io + cr2'^oij i-G- the next- 
to- leading-order terms in Eqs. ( |A6|) - ([A8|) , are obtained by neglecting the L and NL terms. 
They satisfy the equations 
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^ = S,{v, A) , (A22) 

The unique solution of this system which vanishes when A — >^ — oo reads 

01 (t;, A) = dX' dX" S^{v, A") = dX' (A - A') S^{v, A') , (A25) 

J—oo J — oo J — oo 

Pi{v, A) = aiv) 0i(t;, A) + dX' Spiv, A') , (A26) 

J—oo 

mi{v,X) = -a{v)(l)i{v,X) + dX'Smiv,X')- (A27) 



The form of this solution defines the action G*S = J dX! G{X, A') 5'(A') of the retarded matrix 
Green function G(A, A') on any source term ^(A') = (^^(A'), Sp{X'), Sm{X')). As a check on 
the consistency of the iteration method, note that, if one had not imposed any boundary 
conditions, the generic solution ipi would have contained an arbitrary homogeneous solution 
V'hom ("zero mode"). However, it is easy to check that the addition of such a generic i'hom 
would be absorbable in suitable redefinitions of the "seed" functions a{v),Cif,{v), Cpiy), 
Cm{v) entering the background solution Eqs. ( |A3| )-( [A3D . Finally, the explicit form of the 
next-to- leading order solution (with ai = cr2 = 1) is 

^'<-^)^(^^'"'^"^-¥^^"-^-^^ 

A(^;,A) =a0i(t;,A), (A29) 

rm{v, A) = -a Mv, A) + ^, „ e("'+^)\ (A30) 

Note that the coefficients of the terms proportional to e^'^ and A e^'^ contain derivatives with 
respect to v, hence they determine the inhomogeneous expansion near the singularity, while 
the terms proportional to e'^"^^^) ^ are present also in the purely homogeneous case. 

Inserting the double power series expansion ( [A21| ) into ( |A19| ) or ( |A20| ) yields, for each 
contribution ipnm, 
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where the R.H.S. Snm plays, at each order, the role of an effective source term which is 
(in principle) known, being some algebraico-differential expression in the lower-order terms 
ipn'm'i with n' + m' < n + m. Let us make the inductive assumption that Snm is a sum of 
terms of the form 

Sn„^^ = e"("'+^)^e2-^Q,(A), (A32) 

with / < m, where Qi{X) is a polynomial in A of order /. [Here, n,m,l are all integers]. 
It is straightforward to check, from the explicit action ( |A25| )- (|A27| ) of the Green function 
on source terms, that the solution ipnmi corresponding to any term Snmi, Eq- ( |A32|) , in the 
R.H.S. of Eq. ( |A31| ) is again of the form 



i^nml = G * Sr^ml = 6^ ^'^^ ^ 6^^ ^ Q i{X) , (A33) 

where Qi{X) is another polynomial in A, of the same order I. Examining the explicit form of 
the L and NL terms on the R.H.S. of Eqs. (|X9|), (|A13|) and ( A15 ), one can then prove that 



the induction hypothesis consistently propagates at the next order in n + m (still with the 
constraint I < m). 

In conclusion, coming back to the variable r, we find that the inhomogeneous expansion 
near the singularity contains two intertwined series: an expansion in powers of r^"^"*"^) and 
one in powers of mixed with powers of logr. The solution to all orders will be of the form: 

^= E awr"("'+i)r=^'"Pz(logr), (A34) 

n, m, l<m 

where anmi, and the coefficients of the polynomial Pi, are functions of v that in principle 
can be iteratively evaluated. As shown in Sec. |VD| , the full subseries J2n '^noo r'^'^"^"''^^ can be 
exactly summed in terms of the variable x {x 0^) such that = l"^ x^^^ (1 — xY"^ with 
b = (1— a;^)/(l + a^) (see Eqs. ( p.37| ), ( p.38| )). On the other hand, the terms with m 7^ are 
more complicated, as they contain more and more f-derivatives of a{v), C^{v) and Cf3{v) as 
m increases. We note, in passing, that this monotonous increase of f-derivatives shows that 
the series in m (and /) (inhomogeneous expansion) cannot be convergent for generic (non 
analytic) seed functions a{v), C^{v) and Cpiy). 

To complete the proof that the algorithm defining ip solves all the initial field equa- 
tions ( |3.14|) - (|3.17| ), we finally need to check that the m(t>, r) constructed above satisfies the 
constraint ( p.l6| ) involving dm/dv. Let us denote the quantity which must vanish as 
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£ = 



m 



dm" 
dv 



— e 



. dv . 



+ e'' 1- 



2m 



dv 



(A35) 



It is easily checked that, when the other field equations ( p.l4| ), (|3.15|) and (|3.17|) are satisfied, 
8 fulfills the following identity (which could also be derived from the Bianchi identities) 
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— + A{v,r)8 



0, 



A(t>, r) 



i m 



dr 



(A36) 



Noting that the coefficient A{v,r) is bounded near r = {A{v,r) ~ r"^*^'')), and that the 
relation (|5.26|) ensures the vanishing of £^ as r 0, we conclude from the identity ( |A36| ) 
(viewed as an ODE in r with data given at r = 0) that 8 vanishes everywhere. 
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